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FOREWORD 


This report describes work performed by the Sikorsky Division of 
the United Technologies Corporation under Contract NASl-16058 for 
the Structures Laboratory of the U.S. Army Research and Technology 
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The Technical Monitor was Mr. R.L. Tomaine. 
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SUMMARY 


This report describes a Coupled Ro to r/Air frame Vibration Analysis 
which was developed as a design tool for predicting helicopter 
vibrations, and a research tool to quantify the effects of struc- 
tural properties, aerodynamic interactions and vibration reduction 
devices on vehicle vibration levels. The analysis consists of a 
Base Program utilizing an impedance matching technique to repre- 
sent the coupled rotor/airframe dynamics of the system supported 
by inputs from several External Programs supplying sophisticated 
rotor and airframe aerodynamic and structural dynamic representa- 
tion, The Base Program has a stand-alone capability to provide 
assessment of the effects of changes in airframe dynamics and 
vibration reduction devices on vehicle vibration levels. The Base 
Program is not restricted in input support from any particular 
programs and can be executed interactively for many pertinent 
applications. The theoretical background, computer program 
capabilities and limited correlation results are presented in this 
report. Correlation results using scale model wind tunnel results 
show that the analysis can adequately predict trends of vibration 
variations with airspeed and higher harmonic control effects. 
Predictions of absolute values of vibration levels were found to 
be very sensitive to modal characteristics and results were not 
representative of measured values . 
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INTRODUCTION 


The reduction of helicopter vibrations has become increasingly 
important in recent years. Motivations are improvements in 
reliability and maintainability, passenger and crew comfort, 
visual acuity and increased equipment protection. The need to 
reduce vibration levels has influenced the improvement of vi- 
bration analysis methodology. In order to provide an analysis 
capable of accurate vibration prediction, the physical, mathe- 
matical and software attributes listed in Table 1 must be in- 
cluded. A review of the literature (see References 1 to 3) showed 
that published analyses do not possess a complete set of these 
desired attributes. 

From a physical standpoint, the adequate representation of the 
vibration model requires an adequate simulation of the sources of 
vibration, an elastic fuselage, and the elements used to control 
the vibration. Sources of vibration are the air forces acting on 
the rotor, fuselage, and fixed surfaces of the helicopter. 
Vibration control methods include reduction of vibratory loads at 
the sources, modification of response dynamics, and isolation of 
vibration. Limitations occur in one or more of these categories 
and a complete physical model does not appear to be available with 
the scope desired for design studies. 

Regarding the mathematical methodology, the majority of methods 
are limited in their flexibility to represent configurations and 
in their expandability to new configurations since they are based 
on hand-derived equations for the coupled dynamical system, in 
contrast to automated coupling. 

From a software standpoint, there are deficiencies in organiza- 
tion, and internal documentation. Comprehensive modeling capa- 
bility and efficiency require good software to facilitate the 
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TABLE 1 


ATTRIBUTES OF TARGET VIBRATION ANALYSIS 


* Physical 

Comprehensive rotor model 
. Articulated, hingeless, teetering 
. Fixed system 
. Modal 
i . Stick 

Hybrid 

! . Drive Train 

Test Data 

Vibration Control Devices 

Absorbers (hub, blade, airframe) 

Isolators 

Higher harmonic control 
Other vibration sources 
Tail rotor 
! , Empennage 

j . Aerodynamic 

I , 

Rotor induced flow fxeld 
Stall 

Flight modes 

. Hover, forward flight, transition 

i 

!* Mathematical 

General Configuration modeling ability 
j . Steady state forced response, eigensolution, time history 

I . Parameter variation 

* Software 

Efficiency 

User-oriented features - input, plotting 
. Clear Architecture 

. Code clarity 
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addition of new components, and the latitude to exercise and 
improve solution modes and algorithms. The majority of programs 
reviewed were designed before it was recognized that the accom- 
plishment of good software is important for flexibility and ease 
of growth. 

In order to provide the desired attributes a substructure assembly 
method was utilized to simultaneously provide an efficient 
computer code while retaining the required high level of aero- 
dynamic and dynamical representations. The substructure assembly 
method incorporates a stand-alone Base Program representing the 
coupled rotor/airframe system supported by External Programs which 
provide the necessary aerodynamic sophistication. 

This report describes the basis of the coupled rotor/airframe 
analysis, the capabilities and limitations of the analysis, and 
theoretical and practical applications . 
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LIST OF SYMBOLS 


Matrix in complex eigensolution equation (see 
equation (36)), 2N x 2N 

Inverse of system impedance matrix (see equation 
(86) ) , 2N X 2N 

Matrix in complex eigensolution equation (see 
equation (37)), 2N x 2N 

Matrix in Nevm\ark-Beta method (see equation 
(53 ) ) , N X N 

Cosine 

Independent coordinate system damping matrix, 

N X N 

Dependent coordinate system damping matrix, Nj^ x 

Damping matrix for i-th substructure 

Matrix in Newmark-Beta method (see equation 
(48) ) , N X N 

Rotation matrix of direction cosines (see 
Figure 4), 3x3 

Vector of external forces for i-th sub- 
structure, representing forces and moments. 




Components of force applied to a substructure 
connection node in and directions, 
respectively (see figure (2)), N (lb). 


F 


Independent coordinate system force vector, 
N X 1 


F' 


Matrix in Newmark-Beta method (see equation 
(54) ) , N X N 


D 


Dependent coordinate system force vector, x 1 


g 

h 


Independent coordinate system force vector at 

time t in Newmark-Beta method, N x 1 
n 

Acceleration of gravity, m/s^ (ft/s^) 

Time interval in Newmark-Beta method (see equa- 
tion (43 ) ) , s 




Matrix relating rotor hub reaction forces to HHC 
control perturbations (see equation (80)), 12 x 6 


^IIH' ^12H' 
^21H' ^22H 


Submatrices of matrix H (see equation (81)), 
N X 3 


H Independent coordinate system matrix introducing 

effects of HHC control perturbations on system 
vibrations (see equations (84) and (87)), 2N x 6 


^11' ^12' Submatrices of matrix H (see equation (89)), 

^21' ^22 N X 3 
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^IID' ^12D' Dependent coordinate system siibma trices intro- 

^22D ducing effects of HHC control perturbations on 

dependent coordinate system reaction loads (see 
equations (82) and (83)), x 3 

I Identity matrix 

j Performance index for vibration minimization (see 

equation (90)) 

K Independent coordinate system stiffness matrix, N 

X N 




Dependent coordinate system stiffness matrix, N. 


D 


X N 


D 


[k] 


(i) 


Stiffness matrix for i-th substructure 


M 


Independent coordinate system mass matrix, 
N X N 


«D 


m 


m 


m 


xl 

yi 

zl 


[m] 


(i) 


n 


Number of normal modes in a substructure modal 
rep re s ent a t i on 

Dependent coordinate system mass matrix, x 

Components of moment applied to a substructure 
connection node in x^, y^, and z^ directions 
respectively, N-m (Ib-ft) 

Mass matrix for the i-the substructure 

integer step numbers in time-history solution 
(see equation (43)) 
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Index denoting n-th bi filar mass 


N 


N. 


D 


P 


q 

Q 


{r} 


(i) 


r 

r 

r 


xl 

yi 

zl 


^mxl 

^myl 

^mzl 


Total number of independent coordinates in the 
assembled system 

Number of bifilar arms 

Number or rotor blades 

Total number of coordinates in the dependent 
coordinate vector (obtained by summing the 
coordinates for all substructures) 

Matrix in Newmark-Beta method (see equation 

(49) ), N X N 

Modal coordinate 

Matrix in Newmark-Beta method (see equation 

(50) ), N X N 

Vector of reaction forces for i-th sub- 
structure, representing force and moment 
reactions 

Components of reaction force applied to a sub 
substructure connection node in x^, y^, and z^ 
directions, respectively (see Figure (2)), 

N (lb). 

Components of reaction moment applied to a 
siibstructure connection node in x^, y^, and z^ 
directions, respectively (see Figure (2)), 

N-m (Ib-ft) 
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R 

Matrix in Newmark-Beta method (see equation 
(51)), N X N 


Dependent coordinate system reaction load vector, 
1 

s 

Sine 


Scalar coefficient of n-th harmonic of stress or 
pushrod load on a rotor blade, N/m^ (Ib/in^) 


Scalar coefficient of n-th harmonic of stress or 
push rod load on a rotor blade for no hub motion 
and no HHC control input (see equation (94)), 
N/m^ ( Ib/in^ ) 


Matrix introducing the effects of hub displace- 
ment on n-th harmonic coefficient of blade 
stress or push rod load (see equation (94)), 

1 X 12 


Matrix introducing the effects of HHC perturba- 
tions on n-th harmonic coefficient of stress or 
push rod load (see equation (94)) 

t 

n 

time at n-th integration step (see equation 
(43)), s 

T 

Transfer matrix relating HHC inputs to vibratory 
response of the system, (see equations (84) and 
(87)), 2N X 6 

u 

x^ displacement of a connection node (see Figure 
(2)), m (ft) 
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Components of airflow relative to a point on a 
rotor blade, parallel and perpendicular to the 
rotor rotation plane, respectively, m/s (ft/s) 

V displacement of a connection node (see Figure 

2), m (ft) 

w displacement of a connection node (see Figure 

2), m (ft) 


w 


1 


w 


z 


W 


e 


{X} 


(i) 


X 


X 


c 


X 


CO 


Component of airflow normal to the median surface 
of the rotor blade excluding the rotor induced 
flow (see equation (102)), m/s (ft/s) 

Response vector weighting matrix in expression 
for vibration reduction performance index (see 
equation (90)), x 

HHC control angle weighting matrix in scalar 
performance index expression (see equation (90)), 
3x3 

Vector of coordinates for i-th substructure 

Vector of independent coordinates for the as- 
sembled system (=X^), N x 1 

Coefficient vector of cosine harmonic term in 
harmonic expansion of X^ vector (see equation 
(24)), N X 1 

Coefficient vector of cosine harmonic term in 
harmonic expansion of X^ vector for no HHC 
control input (see equations (85) and (86)), 

N X 1 
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X. 


X 


II 


X 


12 


(X^ ) 

' I 'n 


Vector of coordinates for the dependent coord- 
inate system, x 1 

Coefficient vector of cosine harmonic term in 
harmonic expansion of Xj^ vector, Nj^ x 1 

Coefficient vector of sine harmonic term in 
harmonic expansion of X^^ vector, Nj^ x 1 

Vector of independent coordinates for the system 
(= X) , N X 1 

Vector of independent coordinates for system 
after transformation (see equation (14)) 

Vector of independent coordinates for the system 
after transformation (=X^=X) (see equation 

(15)), N X 1 

Vector of independent coordinates at n-th time 
step in time-history solution 


X 


o 


X 


s 


X 


so 


Vector of dependent coordinates resolved to 
global reference axis, x 1 

Coefficient vector of first sine harmonic term in 
harmonic expansion of Xj vector (see equation 
(24) ) , N X 1 

Coefficient vector of sine harmonic term in 
harmonic expansion of X^ vector for no HHC 
control input (see equations (84) and (85)), 

N X 1 
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V 


Vector of coordinates used to obtain the complex 
eigensolution system response (see equation 
(35) ), 2N X 1 



Vector of response variables in the expression 
for the vibration reduction performance index 
(see equations (90) and (91)) 



Impedance matrix relating rotor hub loads to hub 
displacements (see equation (60)), 12 x 12 

^IIH' 

^21H' 

^12H' 

^22H 

Submatrices of matrix z„ (see equation (61)), 

Jq 

6x6 

^IID' 

^21D' 

^12D" 

^22D 

Submatrices of the system impedance matrix 
referred to dependent coordinates (see equation 
(64)), X 



Transformation matrix relating dependent coordi- 
nates to independent coordinates (see equation 
(4)), X N 

Pc 


Coefficient of cosine term in harmonic expansion 
of bifilar mass flap displacement (see equation 
(77)), rad 

Pn 


Flap displacement of n-th bifilar mass (see 
equation (77)), rad 

Ps 


Coefficient of sine term in harmonic expansion of 
bifilar mass flap displacement (see equation 
(77)), rad 

12 





Pi 

P2 

P* 

[V] 


n 


6 


Transformation matrix relating dependent coordi- 
nates Xj^ resolved to substructure local orienta- 
tion axis to dependent coordinates referred to 
global orientation axis (see equation (12))/ 

Nd * Np 

Bifilar mass flap collective displacement, rad 

Transformation matrix relating coordinate vector 
X^ to coordinate vector X^^^ (see equation (14)) 

Transformation relating coordinate vector X^^ to 
coordinate vector X^^ (see equation (15)) 

Newmark-Beta factor 

Modal matrix relating displacements at a node to 
vector of normal mode coordinates (see equation 
(17))/ 6 X M 

Coefficient of cosine term in harmonic expansion 
of bifilar mass lead displacement (see equation 
(74))/ rad 

Lead displacement of n-th bifilar mass (see 
equation (74)), rad 

Coefficient of sine term in harmonic expansion of 
bifilar mass lead displacement (see equation 
(74))/ rad 

Fixed absorber mass displacement, m (ft) 
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0 


Rotor blade collective pitch angle (scalar), 
rad 


Euler angle relating local and global axes (see 
figure (3)) 

e Vector of cosine coefficients of HHC control 

c 

angles (see equation (78)), 3 x 1 

e Vector of sine coefficients of HHC control angles 

s 

(see equation (79), 3x1 




Coefficient of N-th cosine harmonic in blade HHC 
control angle expansion 


(6 ) Coefficient of N~th sine harmonic in blade HHC 

control angle expansion 

Of, 0^, Angular displacements of a connection node (see 

Figure 2 ) 


e, (f) , Euler angles relating local and global axes (see 

Figure (3)) 


Frequency of a normal mode, rad/sec 


Load excitation frequency, rad/s 


Q 


Rotor speed, rad/s 
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Subscripts 


c Coefficient of cosine term in a harmonic expan- 

sion of a variable 

n n-th harmonic of stress or pushrod load 

n-th bifilar mass 
time index 

s coefficient of sine term in a harmonic expansion 

of a variable 

H hub 

Superscripts 

i i-th substructure 

I imaginary part 

n index of last substructure in a set of 

substructures 

R real part 

T transpose 

-1 inverse 

. first derivative with respect to time 
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second derivative with respect to time 
amplitude 

optimum value of a control angle 
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OVERVIEW OF THE COUPLED ROTOR/AIRFRAME VIBRATION ANALYSIS 


The Coupled Rotor/Airframe Vibration Analysis consists of a Base 
Program which is used in conjunction with a set of External 
Programs to achieve a vibration prediction capability for heli- 
copters. A substructure assembly method is used in the Base 
Program to assemble a dynamical model of the helicopter from 
physical helicopter components contained in the Base Program and 
External Program models of the rotor system. The Base Program 
solution modes include steady state forced responses, eigen- 
solutions, time history solutions, and parameter variation capa- 
bility. Rotor data are transferred to the Base Program, which 
functions in a stand-alone mode to facilitate the evaluation of 
parameter variations on vibrations. Stand-alone operation is 
achieved with a small core requirement. The Base Program and the 
External Programs are configured to provide the user with latitude 
to select applications tailored to user interest, aptitude, and 
convenience. If rotor data are on file, the use of the External 
Programs is not required- The Base Program may be used in a 
stand-alone mode also for cases not requiring rotor data - e.g., 
investigation of transraission/isolator/ fuselage response. Input 
to the Base Program is designed to provide an economical format 
requiring input only for selected components. Input to the Base 
Program is not restricted to the External Programs discussed in 
this report. Programs with similar capabilities may be utilized 
provided that input formats are followed. 

Base Program 

The Base Program calculates the mass, damping, and stiffness 
matrices, and forcing function vectors of substructures resident 
in the Base Program code. Input data consist of 3-character 
literal strings which define the components to be included in the 
configuration. These 3-character literal strings trigger input 
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requirements for the components and the forming of the properties 
of the requested components. Input data also define component 
element numbers and points of connection with other components. 
The Base Program uses these data to calculate a transformation 
matrix which synthesizes the coupled equations of motion for the 
assembly from the separate equations for the components. The 
transformation matrix relates the set of coordinates for the 
components of the requested system to a reduced coordinate set 
applicable to the coupled equations of motion. Redundant co- 
ordinates are removed from the coordinates for the requested 
components (the dependent set) to create the coordinates for the 
coupled system (the independent set). Rotation transformations 
relating component local axes to assembled system global axes are 
included in the transformation matrix. The user may elect to 
solve the final equation of motion by one of several methods. 
Component responses are recovered from the solution by multi- 
plication of the solution vector by the system transformation 
matrix. 

The External Programs yield the dynamic and aerodynamic mass, 
damping, and stiffness matrices, and forcing functions of the 
rotor system. These properties are assembled with the properties 
of the components resident in the Base Program by using the 
transformation method described above. 

Figure 1 illustrates the data transfer between the Base Program 
and the External Programs. The arrows indicate the directions of 
the data transfer. Cycling of data between modules is indicated by 
two separate arrows joining program blocks. The character of the 
data transferred is shown in Figure 1. The data transfer from the 
External Programs into the Base Program is accomplished by file 
transfers and the Base Program operates in a stand-alone mode to 
obtain solutions . 
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Figure 1, Coupled Rotor/Airframe Vibration Analysis Data Flows. 



The Four External Programs chosen to support the Base Program to 
form the coupled rotor/airframe analysis are described in a 
subsequent section. 


SUBSTRUCTURE ASSEMBLY METHOD 


The configuration of the Base Program and the data transfer 
between the Base Program and the External Programs are based on a 
substructure approach. The substructure approach synthesizes the 
equations of motion from dynamical models of parts of the system. 
The substructure approach is the key to a flexible modeling 
capability permitting wide latitude in design studies. The 
approach achieves an accurate representation of rotor/airframe 
systems in the framework of an efficient computer program. 


The substructure method employed is the Hurty method of Reference 
4. The coordinates or degrees of freedom of a substructure (also 
called physical component) are physical and generalized displace^ 
ments, such as modal amplitudes. The equation of motion for each 
substructure is expressed in mass, damping, and stiffness matrix, 
and force vector form. Character and treatment of exceptions to 
these forms occurring in the Base Program and the External Pro- 
grams are described subsequently. Properties of the i-th sub- 
structure are [m]^^\ [k]^^\ {f}^^^/ and [r] These 
are respectively mass, damping, and stiffness matrices, and 
external force and reaction force vectors . The s\ibmatrices for 
the substructures are collected into a partitioned diagonal matrix 
equation which represents the system. This partitioned diagonal 
matrix equation is 




( 1 ) 


The matrix M^^ illustrates the typical form of the diagonal 
matrices as : 
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( 2 ) 


”d = 


[m] 


( 1 ) 

[m] 


( 2 ) 


[m] 


(n) 


The vector of coordinates Xj^ illustrates the form of the vectors 
and in equation ( 1 ) . 



The subscript D signifies that equation (1) applies to a set of 
coordinates which may contain coordinate elements which are 
equal to or depend on each other. In other words, Xj^ is a 
"dependent" coordinate set, which can contain redundant elements. 
The coordinate vector X^ contains unconnected coordinates and 
connection node coordinates. The displacement of the mass of a 
simple absorber would be an unconnected coordinate. The dis- 
placement of the base of the absorber would be a connection node 
coordinate because the base of the absorber may be attached to 
another component. When components are connected to each other 
redundant coordinates occur in X^^. 

A force vector, { f } represents external forces and moments 
applied to the i-th substructure. The forces may act at points of 
the structure which are unconnected, like an absorber mass, or at 
a connection node, like the absorber base. Reaction forces act on 
a substructure at points of connection betwen it and other sub- 
structures , and are induced by the other substructures . The 
reaction force vector, does not have any elements asso- 
ciated with unconnected coordinates. Figure 2 shows the sense of 
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Figure 2. 


Displacements and Forces Acting on a 
Connection Node of a Substructure. 
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forces and moments applied to a connection node. Appendix A 
illustrates the terminology introduced for a simple system. 

The synthesis of the equations of motion for the coupled system is 
accomplished by a mapping or transformation relating the dependent 
coordinates, Xj^, to a reduced coordinate set, . Redundant 
coordinates are eliminated by requiring component displacements to 
be equal at connections . The reduced coordinate set may include 
modal coordinates. The transfoimation matrix relating X^^ and X^ 
is denoted by p . The mapping from one coordinate set to the other 
is : 


= P Xj (4) 

Appendix A shows a sample set of coordinates, X^^ and X^ , and a 
sample mapping. 

Substitution of the transformation equation (4) in equation (1) 

. T 

and pre-multiplication of equation <1) by p yields the final 
equation of motion for the coupled system- This is 


MXj + CX 

j + KXj = F 

(5) 

The matrices and vector 

in equation ( 5 ) are 


M = 

pXp 

(6) 

C = 


(7) 

K = 

p%p 

(8) 

F = 


(9) 
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Matrix p is assumed to be a constant function of time, since its 
elements relate displacements at interfaces which are rigidly 
joined, or coordinates which are identical - e.g. internal co- 
ordinates . The condition that reaction forces at connection nodes 
are equal and opposite was used to obtain equation (5). This 
condition reduces to 


= O (10) 

Equation (10) satisfies the requirement that work done by internal 
reactions of the assembled structure in a virtual displacement, 
consistent with the constraints, is zero. 

Equations (6) through (9) are used in the Base Program to synthe- 
size the equations of motion for the coupled system. Appendix A 
illustrates the synthesis of coupled system equations. 

The component displacements may be recovered from the solution 
to equation (5), , by mapping equation (4). Component veloc- 
ities and accelerations are derived similarly. With and its 
derivatives known, connection node or interface reactions from 
equation (1) may be expressed as 


Rotor hub shears and moments are sample interface reactions. 
Appendix A illustrates the recovery of interface reactions. 

Composition of Transformation p 

The transformation matrix p is expressed as the product of three 
transformations which assemble coupled systems from components 
whose properties are defined in component local axes, and allow 
for the use of modal coordinates. These transformations are 
described in the following subsections. 
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Transformation from Local to Global Axes, 

The transformation from dependent coordinates resolved to local 
axes with arbitrary angular orientations to dependent coordinates 
resolved to a global reference axis is 


^ = Po^o (12) 

Dependent coordinate sets and can include coordinates which 
are redundant. Redundant coordinates occur when components are 
connected to each other at connection nodes . Equation ( 12 ) 
resolves coordinates at connection nodes to the same system and is 
a prerequisite to the elimination of redundant coordinates from 
X^. Figure 3 illustrates the angular orientation of a local axis 
system relative to a global system. The local system is obtained 
by three successive rotation displacements of axes through Euler 
angles 0 , (j> , and . 


The rotation transformation relating linear displacements in the 
local and global axes is 

/ u \ 

= [E] < V > (13) 

Local ^ Global 

Matrix [E] is a 3 X 3 matrix of direction cosines each of which is 
a function of Euler angles 6 , 41 , 4 / . The same transformations 

apply also when small angular displacements 6 ^, 62 / 63 are sub- 

stituted for u, V, w. Linear and angular velocities and ac- 
celerations are mapped by equivalents of equation (13). Figure 4 
defines the elements of matrix [E] . Appendix A shows the ap- 
plication of a rotation transformation to a simple absorber. The 
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absorber mass displacement, 6, is an unconnected coordinate and is 
retained as a coordinate in the final system. It is not necessary 
to resolve 6 to a global system because unconnected coordinates 
do not have to be eliminated by equating displacements at connec- 
tions. Unconnected coordinates are assumed to be invariant to 
rotation transformations from local to global axes. 

Transformation from Dependent Coordinates to Physical Domain 

Independent Coordinates, 

The transformation from dependent coordinates referred to a global 
reference axis to coordinates from which redundant coordinates at 
connection nodes have been removed is 

Xq = Pi Xj^ (14) 

Coordinate set X^^, is referred to as the physical domain inde- 
pendent coordinate set since connection node displacements are 
physical displacements. The determination of is based on a 

convention in the computer program specifying the order of co- 
ordinates in the dependent sets X and . Coordinates in X are 
ordered according to ascending component element numbers assigned 
by the user in the input specification to the Base Program. 
Appendix A illustrates this organization in X^. In coordinate set 
Xj^ , unconnected coordinates are listed first in order of ascend- 
ing element number. Connected coordinates are listed after 
unconnected coordinates in order of ascending connection mode 
numbers . Redundant coordinates are eliminated by first retaining 
coordinates for a component with the lowest element number, and 
then eliminating redundant connected coordinates with higher 
element numbers from a connection node. Appendix A shows the 
ordering of coordinates in X^^^ . The logic yielding the elements 

of p^ is based on the ordering of coordinates in X and X^- 
1 oil 
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described above. Appendix A shows the transformation matrix 
corresponding to the illustrations given for and . 

Transformation from Physical Domain Independent Coordinates to 
Modal Domain Independent Coordinates, 

The transformation from physical domain independent coordinates, 
X T , to coordinates which include modal coordinates, X^^ , is 

^11 ~ ^2 ^12 ( 15 ) 

This transformation is required when modal structures are used. 
Since Xj 2 1® final coordinate set for the system to be solved 

it is symbolized also by X^. , i.e. 

Xj = ( 16 ) 


The coordinate set X^^ Is obtained from X^^ by expressing 
connection node coordinates in terms of modal coordinates. For 
example, rotor/airframe hub coordinates may be expressed in terms 
of airframe modal coordinates. Coordinate set X^^ contain 
connection coordinates which have not been expressed as modal 
coordinates, and will in general include unconnected coordinates 
for the substructures. Coordinate set X^^ is referred to as the 
independent coordinate set in the modal domain. 

Connection node displacements are expressed in terms of modal 
coordinates {q}. The linear and angular displacements at a con- 
nection node referred to a local axis for the modal structure are 
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u 


V 

w 

e. 


= [y] {q} 


Local 


(17) 


In equation (17), [y] is a 6XM modal matrix. The vector of normal 
modes of the substructure has dimensions MXl, i.e., the sub- 
structure equations of motion have M normal modes. Small angular 
displacements are assumed. Velocities and accelerations are 
mapped by the equivalents of equation (17). The local axis system 
displacements for the modal structure are related to the global 
axis displacements for the final system by 



Local 


Global 


(18) 


In equation (18), [E] is the matrix of direction cosines introduc- 

ed in equation (13) with elements defined in Figure 4. To obtain 
the connection node displacements referred to the global system in 
terms of modal coordinates, the inverse of the rotation transforma- 
tion [E] is expressed as: 


-1 T 

[E] = [E]^ 


( 19 ) 
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Global 




:ation Transformation Matrix [E]. 


I 



The solution to eguations (17) and (18) is 



[E] I O 


O 


I T 


[vHq] (20) 


Combining equation (20) with equation (15) yields the matrix 
relating the coordinates in to modal coordinates and other 

coordinates in Equation (20) is used at each connection node 

where physical displacements are to be expressed in terms of modal 
displacements. In general, the matrices [E] and will differ 

among the nodes because different modal substructures may be 
joined to different nodes of the system, and because of different 
component orientations. Appendix A gives examples of 

^ 2 * 


Composite Transformation p Relating Dependent Coordinates to 
Independent Coordinates 

From equations (12), (14), (15), and (16) it is found that 


^ = Po Pi P2 


( 21 ) 


Comparison of equations (4) and (21) yields 


P 


P-, P 


( 22 ) 
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This transformation relates the dependent coordinate set, for 
the assembly of substructures to the final independent (non- 
redundant) coordinate set, X^, used to solve the equations of 
motion. Equation (22) is used to construct the p -transformation 
in the computer program. Appendix A gives an example of the 
composite transformation. 


SOLUTION MODES 

Four solution modes of equation (5) can be obtained. Each method 
has different attributes affecting information provided, speed of 
execution, and computer program core requirements. These solution 
modes are: 1) steady state forced response, FRl, 2) eigensolution 
mode for systems without damping, EGl, 3) eigensolution mode for 
systems with damping, EG2, and 4) time history mode yielding 
system response as a function of time, THl . 

In addition to the solution modes described above, two more opera- 
tional modes of the program may be invoked. These modes are 1) 
parameter variation mode, PVl, designed to provide a convenient 
selection and variation of parameters in design studies, and 2) 
general program control mode, GEN, presently configured to control 
display of program output and graphics options . 

Table 2 lists the character names associated with each solution 
mode and operational mode available. 


A description of each solution and operational mode is given 
below: 
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TABLE 

2 - BASE PROGRAM SOLUTION MODES AND OPERATIONAL MODES 

■ 

MODE 

DESCRIPTION OF MODE 

FLAG 


FRl 

FORCED RESPONSE 1 - Forced response solution mode 
yielding harmonics of steady state response 

EGl 

EIGENSOLUTION 1 - Eigensolution mode for systems 
without damping 

EG2 

EIGENSOLUTION 2 - Eigensolution mode for systems 
with damping 

THl 

TIME HISTORY 1 - Time history solution mode using 
Newmark-Beta method 

PVl 

PARAMETER VARIATION 1 - Operational mode facilitat- 
ing variation of parameters for design studies 

GEN 

GENERAL - General control mode providing general 
control for a case 




steady State Forced Response (FRl) 


The harmonic steady-state forced response mode FRl, yields 
harmonic coefficients of response to a steady harmonic excit- 
ation. This solution mode is useful for design studies by virtue 
of execution time, small system core requirements, and importance 
of steady state vibrations. 

The forcing function which appears in equation (5) is assumed to 
be 


F = F cos u)j^t + F sin U3^t (23) 

c f s f ' ^ 

In equation (23), u)^ is the forcing frequency and F^ and F^ are 

N X 1 vectors, where N is the dimension of vectors in equation 
(5). The harmonic steady state forced response is also assumed 
periodic . 


= X cos uj^t + X sin uj^t (24) 
I c f s f ' ^ 

Substitution of equations (23) and (24) into (5) and equating the 
coefficient expressions of cos m^t and sin lu^t to zero yields the 
2N X 2N algebraic equation 
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( 25 ) 


■U3^2m+K j U)^C 


■u)^C 


) -U)^2m+K 


The solution to equation (25) is 



F 


c 


F 


s 



-U)^2m+K I 

X I 


to 




-lu^c I -U)^2jyi+K 


1 


F 

c 


F 

s 


(26) 


Determination of vibratory response requires one inversion, 
equation (26). Solution speed is much higher than a time history 
solution and the method is well suited to parametric design 
studies . 


Equations (25) and (26) apply to systems expressible in mass, 
damping, and stiffness form. These matrices may contain aero- 
dynamic contributions such as in the matrices derived from aero- 
elastic analysis E927. In aeroelastic analysis G400, rotor 
characteristics, including aerodynamic effects, are represented as 
a hub force and moment excitation vector and a hub impedance 
matrix. Rotor induced impedance matrices must be added to the 
system matrices in equation (26) to obtain the vibratory res 
ponse. The steady state solution is 


X 


X 




-1 


u)^C+Zi2 


U 4 


— u) ^C+Z^ - 
j_ *"f^ ^21 


I “‘"i M+KtZ22 


(27) 
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Forcing function vectors F and F in equation (27) contain rotor 

o s 

contributions. Impedance matrices in equation (27) apply to the 
assembled system with independent coordinates, X and X . Trans- 
formations similar to those used to transform mass, damping, and 
stiffness matrices in equations (6) through (8) yield the imped- 
ance matrices in equation (27). The subject of mapping rotor 
impedance matrices to independent coordinates is treated in a 
subsequent section. 

The system vibratory response can be expressed in units of accel- 
eration of gravity by 

X^/g = -(u)^2/g)X^ (28) 

X^/g = -(u)^2/g)x^ (29) 

The computer program displays the amplitude and phase of each 

coordinate in X^ . This information is extracted from the 

elements of X and X . The response of each coordinate is assumed 
c s 

to be a cosine function. The standard line-printer output does 

not display X and X . A debug option may be flagged to obtain X 
o s o 

and X^ . Amplitude and phase data only are stored for plotting. 

Eigensolution for Systems Without Damping (EGl) 

The eigensolution mode for systems without damping yields the 
frequencies and mode shapes of structures without damping. 

Matrices must be constant with respect to time. The effects of 
changes in physical properties on mode shapes and frequencies can 
be determined. These effects are used to aid the interpretation 
of the effects of design variables on normal modes and vibratory 
response. This solution mode has a potential capability, not 
implemented in the present version of the base program, to trans- 
form substructure mass, damping, stiffness and force vector 
properties to normal mode coordinates . Normal mode coordinates 
reduce the problem to one requiring fewer coordinates . 
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Equation (5) is specialized to 

MXj + KXj = 0 (30) 

The solution is assumed to be periodic 

Xj. = Xe^“^ (31) 

Substitution of equation (31) into (30) yields the NxN algebraic 
equation 


(-u>2m + K) X = 0 (32) 

Solutions to equation (32) yield the eigenvalues, u)2, and cor- 

responding eigenvectors, X. 

Eigensolutions for Systems with Damping (EG2) 

The eigensolution mode for systems with damping is used to provide 
information on the stability of coupled systems. This solution 
mode is restricted to hover. Coupled systems may include rotor 

representations which embody dynamic and aerodynamic matrices 

derived from rotor aeroelastic analysis E927. Stability informa- 
tion may be used to supplement the steady state vibratory response 
data for a system to provide more data on system attributes. 
Matrices must be constant with respect to time. 

Equation (5) is specialized to 

MXj + CXj + KXj = O (33) 

Reduced to the first order system equation (33) becomes 
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AY - BY = O 

This is accomplished by the substitution 




. \ 




\^l( 


Y 

_ s 

lx ( 



( 

I ) 


equation 

(34) 

are 




1 

M 

A 


M 1 

C 



"m 1 

0 

B 

= 


— 



[o 1 

-K 


(34) 


(35) 


(36) 


(37) 


The eigenvalue problem is obtained from equation (34) by assuming 


Y = Ye (38) 

which yields the 2N X 2N algebraic equation 

(A A - B)Y = O (39) 

The modal frequency in hertz and modal damping in fractions of 
critical damping respectively are 


w = A^/(27t) 


(40) 


with 


4 = -A /(A. 2 + A 2)^ 
^ r^ ' 1 r ' 


A = A + i A . 

r 1 


(41) 


(42) 
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representing a complex eigenvalue. The real and imaginary com- 
ponents of Y are displayed for each mode. A stable mode corres- 
ponds to a positive value of ^ (negative value of 


Time History Solution Mode (THl) 

The time history solution mode THl, yields the response of a 
system as a function of time. Transient responses, steady 
responses, and the stability of systems can be determined from the 
time-history behavior of a system. Responses of nonlinear systems 
can be determined. The time-history mode is able to, or has the 
potential to, supply all the information provided by forced 
response and eigensolution modes, in addition to the analysis of 
transient and non-linear data. A drawback of the method is its 
high execution time in comparison with the other methods making it 
comparatively costly and inefficient for design studies. 

The solution algorithm yielding the time history response is the 
Newmark-Beta method described in Reference 5. The method is a 
finite difference method. Displacement responses are obtained 
from displacements known at prior times, and from data defining 
the magnitudes of known external forces acting on the system. 

The vector of coordinates satisfying equation (5) at time step n 
is denoted The corresponding time is 

t = t - + nh (43) 

n n-1 ' ’ 

In equation (43) h is the step size. The integer step number n 
ranges from one to the number of steps in the calculation. Time 
t is the initial time. Corresponding initial conditions are 

(Xi>o 
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When (Xj)q and (X^)^ are specified the displacement response at t^ 
is found from 

(Xj)^ = D"^ + Q(Xj)^ + p*h2F^ + RF^] (44) 

Matrices D, P, Q and R are functions of M, C, and K in equation 

(5), time step h, and factor p* (see equations (48) through (52)). 

They have dimensions NXN where N is the size of in equation 

(5). The Factor p* is the Newmark-Beta Factor. Values of p* 

ranging from O to 0.25 can be input to the Base Program. The 

factor p* is used to control the variation of acceleration assumed 

in the time interval (see remarks in Reference 5). A numerical 

solution with p* = 0.25 is unconditionally stable. Force vectors 

F-. and F are known external forcing functions at t and t^ . 

1 o o 1 

After the solution known, successive solutions are 

obtained from the recursion formula 

+ P*h2 (F^_^^ + (l/p*-2)F^ + (45) 

The NxN Matrices B and F' are known functions of M, C, and K, step 
size h, and factor p* (see equations (53) and (54)). 


Equation (45) may be used also to restart solutions from solutions 
calculated up to the time of restart. 


The additional assumption is made in the Base Program code that 
the force F^^^ linearly extrapolates the forces F^ and F^_^. For 
equal intervals of time h this assumption reduces to 


n+1 


= 2F„ - 


n-1 


(46) 
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and equation (44) becomes 




Equation (46) allows for the inclusion of aerodynamic forces which 

r ** 

are functions of Xj , and Xj . The present version of the Base 
Program does not model these aerodynamic forces. A technique to 
include these forces is to assume F to be a function of X^ and its 
derivatives at a prior time. Equations (46) and (47) are consis- 


tent with this assumption. 

Matrices in equations (44), (45), and (47) are 

D = M + (h/2)C + p*h2K (48) 

P = M + (h/2)C - (1/2-p*) h2K 

-(l/4-^*)h3CM"^K (49) 

Q = [M - (1/4 - ^*) h2CM"^C] h (50) 

R = [(1/2 - p*)I + (1/4 -p*) hCM“^ ] h2 (51) 

I = unit matrix (52) 

B = 2M - (1-2P*) h2K (53) 

F’ = M - (h/2)C + p*h2K (54) 


The dimensions of matrices in equations (48) through (54) are NXN 
where N is the dimension of the vector of independent coordinates 
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Parameter Variation Mode (PVl) 


The parameter variation mode, PVl, is used to identify and specify 
a range for a parameter which is to be changed in design studies. 
This operational mode provides a convenient and general mode of 
selection of parameters to be varied among the substructures 
available in the Base Program. Up to ten parameters from any of 
the substructures may be varied simultaneously. 

General Control Mode, GEN 

The general control mode, GEN, is designed to provide a basis for 
the application of general controls to a case such as line printer 
output and graphics selection. 


ASSEMBLY OF ROTOR IMPEDANCE MATRIX 

The purpose of this section is to derive the procedure for the 
assembly of a rotor impedance matrix with mass, damping, and 
stiffness matrices and forcing functions of other substructures. 
The mapping of the rotor impedance matrix to independent coordi- 
nates for the coupled system is shown to be the same as for other 
substructure properties. The rotor impedance matrix must be 
partitioned to accomplish the mapping. 

The effects of hub displacement are included in the coupled 
rotor/airframe system by representing the rotor effects as the sum 
of a hub excitation for a fixed hub and a correction involving the 
rotor impedance matrix multiplied by hub displacement . This 
representation enables rotor properties to be calculated prior to 
execution of the Base Program. This formulation maintains the 
desired small core and execution time for the Base Program, 
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To obtain a rotor hiib excitation vector and a rotor hub impedance 

matrix, the rotor is first analyzed for a fixed hub in the aero- 

elastic analysis - e.g. G400. The six hub forces and moments 
acting on the fixed hub are 

{f} jj = {f^} cos m^t + {fg} jj sin u)^t (55) 

Vectors {f } „ and {f } „ are 6X1 component vectors of cosine 

and sine terms in equation (55). For a N-bladed rotor with 
identical blades, frequency is equal to nNQ, where n is an 

integer and Q is the rotor speed. Forces and moments at frequency 
are the only vibratory forces that can be transmitted to the 
non-rotating axis at the hub. 

To calculate elements of the rotor impedance matrix in aeroelastic 
analysis G400, the rotor hub is perturbed harmonically in the 
fore-aft direction by a displacement u^cos u)^t. This displacement 
causes a force perturbation acting on the hub . The reaction force 
acting on the hub can be expressed as 

{r} = (f^l jj cos w^t + {fg} ^ sin lu^t 

-{z^^} u^ cos u)^t - sin lu^t (56) 

The 6X1 reaction force is equal to the constrained hub forces 
augmented by terms which are proportional to the hub fore- aft 
displacement, u^. The 6X1 column vectors and are 

columns of the impedance matrix. Elements in these vectors are 
assumed to be constant. In a similar way, lateral perturbations 
of the hub, v^cos w^t, are introduced to derive impedance columns 
(z _} and {z _} , This procedure is repeated for all possible 

displacements at the hub, including displacements proportional to 
sin u)^t. The complete rotor impedance matrix and the expression 
for the hub reactions including all hub displacements is obtained. 
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The following definitions express harmonic coefficients for the 
hiib reactions as functions of the impedance matrix. The 6X1 hub 
reactions are assumed to be 


{r} jj = {r^} g cos uj^t + H sin u)^t 

The 6X1 hub displacement perturbation vector is 

{x} = {x^} jj cos uj^t + {Xg} ^ sin u>^t 

with 


(57) 


(58) 


{x} ^ = {u, V, w, e^, e^} 


(59) 


Expressions similar to equation (59) define {x } and [x } „ 

C xl S il 

Equating the coefficient expressions for cos u)^t and sin u)^t 
in the reaction force equation to zero yields 





h/ 


T9 

) H 

H 


h( 



H 


(60) 





^IIH 

1 

^12H 




(61) 




^21H 

1 

( 

^22H 





Matrix 

z„ is a 12 X 

ri 

12 rotor hub 

impedance 

matrix 

composed of four 

6X6 

submatrices 

^IIH' 

^12H' 

^21H' ^22H 

as shown 

in equation 

(61). 

The hub reaction 

force 

in 

equation 

(60) 

is 

composed of 


fixed hub forces {f } and {f^} „ 

C xl S M 


and a correction proportional 


to the impedance matrix multiplied by the hub displacement 
vector. Substituting equation (61) into (60) and rearranging the 
terms yields 
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( 62 ) 


IIH 

1 z 

, 12H 

21H 

1 

1^22H 


^*■ 3 ^ h) ( H 


To develop the algoritlm for assembling impedance properties in 
equation (62) with properties of other substructures expressed in 
mass, damping, and stiffness form, a rotor impedance substructure 
is introduced. This has null mass, damping, and stiffness 
matrices, and non-zero hub force and moment excitation and rotor 
hub impedance. The vector of the dependent coordinates is 
expressed as 

External forces Fj^ and reactions are expressed similarly. 
Substructure mass, damping, and stiffness matrices are collected 
with impedances to yield 


- ^ ^IID j 

^f^D ^12D 

_| 


- “fS + ^21D I" 

*^f^^ ^ ^ ^22D 


(^D>s) ((^>s) 



(64) 


Vectors (Xp)^ and (X^^) 
{x } „ from equation 

S xl 

matrices are diagonal . 


contain rotor hub coordinates {x } „ and 
s ‘ c^ H 

(62). In equation (64), the partitioned 
Matrices Mj^, and Kj^ are typified by 


45 



( 1 ) 


= 


[m] 


[m] 


( 2 ) 


[O] 


H 


[m] 


(n) 


(65) 


Zero diagonal submatrices [0]„ are associated with coordinates 

rl 

{x } „ and {x } „ which have no mass, damping, or stiffness 
c n s n 

matrices in the rotor impedance substructure defined by equation 
(62). The impedance matrix 2^^^^ in equation (64) is 




IID 


[O] 


( 2 ) 


'IIH 


[O] 


(n) 


( 66 ) 


The matrix in equation (66) is zero except for a 6 X 6 diagonal 

element Zt-,„ associated with hiib vector {x Element 

IIH c H IIH 

derives from the impedance substructure defined in equation (62). 
The forms of impedance matrices ^21D' ^22D equation (64) 

are similar to equation (66). 


Substructure RES which models a rotor hub excitation force and a 
hub impedance matrix incorporates the properties described above. 
Each 6X6 mass, damping, and stiffness matrix for RES is zero. 
Substructure impedance elements and hub excitation forces are 
transmitted from external aeroelastic analysis G400. 

The derivation of the transformation to independent coordinates is 
accomplished by substituting equation (6S) for X^^ and equation 
(24) for Xj into equation (4). This yields 
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P X 

c 

(67) 


= 

P X 

s 

(68) 


Substitution of equation (67) and (68) in equation (64) and pre- 
multiplication of equation (64) by 


H 

CO. 

1 

O 

L ° i 

T 

P 


(69) 


yields equation (27) for the vibratory response of the coupled 
system given by 


X 


X 




u)fC+Zi2 


-1 






(27) 


Impedance matrix mappings yielding the matrices in equation (27) 
are 


'll 

^ ^IID ^ 

(70) 




12 

^ ^12D ^ 

(71) 




12 

^ ^21D ^ 

(72) 

22 

^%2D P 

(73) 
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In the computer program implementation the zero mass, damping and 
stiffness properties of substructure RE3 are assembled with the 
properties of other substructures to form and in 

equation ( 1 ) . These matrices are mapped to the independent 

coordinate system by equations (6) through (8). The forced 
response solution mode is invoked and matrices in equation (26) 
are formed. Impedance matrices are mapped to independent co- 
ordinates by equations (70) through (73) and are added to the 
matrices of equation (26) to form equation (27). Equation (27) 
yields the vibratory response of systems including a rotor imped- 
ance model . 

Hub reactions can be recovered from equations (67), (68), and (64) 
after the coupled system solution is found from equation (27). 

The RES impedance substructure is consistent only with a steady 
state forced response mode of solution, FRl. The forcing fre- 
quency in FRl must be the same as the forcing frequency used to 
derive rotor impedance properties. This frequency will be an 
integer multiple of rotor speed and blade number (tu^ = nN Q). 

BASE PROGRAM SUBSTRUCTURES 


Substructures available in the Base Program listed in Table 3 are 
described below. Descriptions are provided of substructure 
function, coordinates, properties, and assumptions. 

Substructure FAl 

The function of substructure FAl is to simulate a vibration 
absorber which may be attached to any other substructure in the 
fixed - i.e., non-rotating-system. Substructure FAl has seven 
coordinates. These are absorber mass displacement, and six 
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TABLE 3 - BASE PROGRAM SUBSTRUCTURES 

COMPONENT 

FLAG 

DESCRIPTION OF COMPONENT 

FAl 

FIXED ABSORBER 1 - Single mass, spring and damper 
vibration absorber. 

GFl 

GENERALIZED FORCE 1 - Force component allowing speci- 
fication of harmonic force and moment excitation. 

BMl 

BEAM MODEL 1 - Uniform elastic beam segment. 

CNl 

CONSTRAINT COMPONENT 1 - Component allowing elimina- 
tion of coordinates (degrees of freedom) from a 
connection node. 

RE2 

ROTOR ELASTIC 2 - Model of aeroelastic rotor expressed 
in mass, damping, and stiffness form. Rotor data are 
transmitted to the Base Program by rotor aeroelastic 
analysis E927. Restricted to hover. 

MSI 

MODAL SUBSTRUCTURE 1 - Normal mode representation of a 
dynamical structure. 

BFl 

BIFILAR 1 - Horizontal linear bifilar vibration 
absorber for reducing rotor inplane vibrations. 

BF2 

BIFILAR 2 - Vertical linear bifilar vibration absorber 
for reducing rotor out of plane vibrations. 

ISl 

ISOLATOR 1 - Nodal isolator employing two anti-resonant 
bar weights . 

RE3 

ROTOR ELASTIC 3 - Impedance model of aeroelastic rotor 
in hover or forward flight. A hub force and moment 
excitation vector and a hub impedance matrix are trans- 
mitted to the Base Program by rotor aeroelastic analysis 
G400. Vibratory response of a coupled system is found 
using solution mode FRl . This component has a cap- 
ability to obtain optimal HHC control angles for vibra- 
tion reduction. 
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Figure 5. Fixed Absorber FAX. 
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displacements of the base consisting of three linear base dis- 
placements and three angular base rotations. The mass of the 
absorber is free to move perpendicular to the base and does not 
move parallel to the base. Viscous damping is assumed. Figure 5 
illustrates substructure FAl properties. Initial mass displace- 
ment and velocity can be specified for transient response ana- 
lysis . 


Substructure GFl 

Generalized Force 1, GFl, specifies a harmonic force and moment 
excitation acting at a point on a structure. The GFl substructure 
has six coordinates, comprising three linear displacements and 
three angular rotations at a connection node. Forces and moments 
are resolved to a non-rotating axis system. The user specifies 
the cosine and sine coefficients of these forces and moments to 
define a harmonic excitation of known magnitude. The mass, 
damping, and stiffness matrices for substructure GFl are null 
matrices . Substructure GFl may be used to represent an excitation 
for forced response (FRl) or time-history (THl) solution modes. 
Substructure GFl may be included in eigensolution (EGl and EG2 ) 
input streams without affecting the final eigensolution. Figure 6 
illustrates substructure GFl properties. 

Substructure BMl 

Beam Model 1, BMl, is a uniform elastic beam. This substructure 
has 12 coordinates consisting of three linear displacement and 
three angular rotations at each end of the beam. The beam may be 
connected to other substructures. The uniform mass of the beam, 
m^, is approximated by two concentrated masses equal to ^^72 at 
its ends. The beam is assumed to have zero mass moments of 
inertia associated with bending and torsion rotations. Structural 
damping is assumed to be zero. Figure 7 illustrates substructure 
BMl properties. 
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Figure 6. Generalized Force GFl. 
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Constraint Substructure CNl 


Constraint substructure CNl may be used to eliminate or constrain 

coordinates at a connection node. There are six coordinates at 

each connection node (see Figure 2). The coordinates consist of 

three linear displacements and three angular rotations. The user 

may constrain any five of the six coordinates at a connection 

node. At least one coordinate among the six must be retained. 

The constraint component is a specialized modal substructure (see 

Substructure MSI) whose modal matrix contains ones and zeroes. 

The occurrence of a one in a modal matrix row causes a connection 

node coordinate to be replaced by a modal coordinate which is 

identical to the connection coordinate. The occurrence of zeroes 

in a modal matrix row eliminates a connection node coordinate. 

Equation (20) explains how the connection node coordinates are 

T . . 

eliminated. Matrix [E] , containing the rotation transformation 
from local to global axes is replaced by the identity matrix I. 
The number of coordinates in vector {g} is selected to be equal to 
the retained coordinates at the connection node. The coordinates 
in vector {q} are identical to coordinates to be retained at the 
connection. 

Rotor Aeroelastic Substructure RE2 and External Program E927 

Rotor aeroelastic substructure RE2 is used to represent a multi- 
blade aeroelastic rotor in hover. This component identifies rotor 
matrices transmitted from the rotor aeroelastic analysis E927 and 
arranges the rotor data for assembly with other components. 
Standard application of RE2 is in conjunction with eigensolution 
mode EG2 . 

Table 4 lists RE2 coordinates, expressed as output labels dis- 
played by the Base Program. Use of rotor coordinates in the rotor 
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TABLE 4 

- AEROELASTIC ROTOR RE2 COORDINATE LABELS 

COORDINATE 

DESCRIPTION 

OUTPUT LABEL 


QFT1-QFT4 

Blade bending collective modal 
amplitude^ 

THTR 

Blade pitching^^^ 

THTE 

Blade torsion^^^ 

BETA 

Flapping collective mode^^^ 

GAMA 

Leading collective mode^^^ 

QTIS, QTIC 

Blade bending cyclic sine and cosine 

QT2S, QT2C 

modal amplitudes ^ 

QT3S, QT3C 
QT4S , QT4C 


THRS , THRC 

Blade pitching sine and cosine cyclic 
amplitudes 

THES , THEC 

Blade torsion sine and cosine cyclic 
amplitudes ^ 

BETS , BETC 

Blade flapping sine and cosine cyclic 
amplitudes ^ ^ ^ 

GAMS , GAMC 

Blade leading sine and cosine 
amp 1 i tudes ^ ^ ^ 

X, Y, 2, THTX, 

Hub displacements . See Figure 2 for 

THTY, THTZ 

positive directions 


NOTES : a ) Positive flatwise bending up and edgewise bendin 

forward. 

b) Positive nose down. 

c) Positive up. 

d) Positive forward. 
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equations, and the assvimption of identical blades, equally dis- 
posed on the rotor disc, enables a constant coefficient equation 
set to be derived ( see Reference 6 ) . Rotor coordinate selection 
is made in the input specification to E927. Rotor aeroelastic 
analysis E927 transmits to a file the number of coordinates, flags 
defining coordinates selected, and corresponding mass, damping, 
and stiffness matrices. Matrices embody dynamic and aerodynamic 
contributions . Numerical input data following the RE2 flag 
consist of hub connection node number and hub orientation angles. 
Processing by the Base Program includes the addition of a dummy 
yaw coordinate. Zero mass, damping, and stiffness matrices are 
created for this coordinate. This yields a six coordinate h;ib 
consistent with coordinate representation for other substructures. 
Six hub coordinates enable the rotor component to be coupled to 
other components without changing Base Program assembly logic. A 
valid physical model of rotor yaw motion is not achieved. Null 
rotor matrices are associated with this displacement. The user 
must eliminate the hub yaw coordinate from the system model. This 
may be accomplished with constraint component CNl . A valid model 
would occur also if fuselage normal modes in component MSI have no 
hub yaw displacement. 


Modal Substructure MSI 

Modal Substructure MSI represents a dynamical system employing 
normal modes of a substructure. Each substructure of MSI type is 
limited to the use of one normal mode. The user must supply input 
defining the modal generalized mass, damping ratio, and frequency 
of the normal mode. These properties can be generated from 
analysis or from test data. Processing in the program forms 
damping and stiffness of the mode from the input data. Each 
substructure of type MSI is allowed to have up to five nodes at 
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which node number, structure orientation angles, and elements of 
the modal matrix for the normal mode are defined. At each such 
node, the elements of the modal matrix are three linear displace- 
ments and three angular rotations for a unit value of the normal 
mode amplitude. The directions of the displacements at a node 
must be consistent with the convention shown in Figure 2 . The 
specification of the same node number between two substructures of 
any type leads to the connecting of the two substructures at the 
node. More than one modal substructure of type MSI is connected 
to the same node when it is desired to represent the effects on a 
node of several normal modes. Nodes on MSI substructures may be 
used to monitor airframe responses at up to five points. If a 
physical substructure is not connected to a node at which re- 
sponses are desired, these responses may be displayed by con- 
necting a GFl substructure to the node of interest, with zero 
forces and moments specified. This GFl substructure contributes no 
matrices or forces to the system. It introduces physical dis- 
placements at the node which become members of the set of depend- 
ent system coordinates which are used to monitor the node response . 
Initial modal amplitudes and velocities can be specified to 
analyze system transient response. 

Horizontal Linear Bi filar BFl 

Horizontal linear bifilar BFl represents a bifilar vibration 
absorber designed to reduce N per rev vibrations in the plane of 
rotation of a helicopter rotor. The bifilar absorber is assumed 
to consist of N arms equal in number to the rotor blades. Each 
arm holds a pendular mass free to lead or lag. The arms are 
assumed to be equally disposed and each holds an identical ab- 
sorber mass. The bifilar absorber rotates at the helicopter rotor 
speed Q. Figure 8 is a schematic of one arm of the bifilar. Each 
bifilar mass lead-lag displacement may be expressed as the sum of 
bifilar lead-lag collective and cyclic displacements. The lead- 
lag displacement of the n-th bifilar is 
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N MASSES. ONLY ONE MASS SHOWN 

y =(y COS<l/ + SIN )(2/N) 
c s 


Figure 8. Horizontal Linear Bifilar BFl. 
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(74) 


= Vq/N + (Y^cosij;^ + YgSin 4»nH2/N) 

= (27t /N) (n-1) + (75) 

i{; = Qt (76) 

Equation (74) introduces the bifilar collective and cyclic coordin- 
ates in the equations of motion of the bifilar. The angle is 

the azimuth angle of the n-th bifilar arm. The equations of 
motion for each bifilar consist of N equations defining the 
equilibrium of moments about the hinge point of each pendular 
mass. These N equations are supplemented by six equations defin- 
ing the equilibrium of the bifilar hub in three translational and 
three rotational directions. When equation (74) is substituted 
into these equations, an equation set with eight coordinates is 
obtained. The assumption is made that = 0- The eight co- 

ordinates consist of two bifilar mass sine (Yg) ^nd cosine 
cyclic amplitudes and six hub displacements shown in Figure 8. 
The reduced equation set has constant coefficients. Restrictions 
on the configurations that can be analyzed are bifilars which have 
three or more equally disposed and equal masses. Pendulum dis- 
placements are assumed to be small enough to permit linearizations 
of the equations. Solutions are obtained in terms of the cyclic 
amplitudes y and y , as functions of time. Solution modes that 
may be used with BFl are FRl, EGl, EG2, and THl . Initial values 
of bifilar mass cyclic mode amplitudes can be specified for tran- 
sient analysis. 

Horizontal bifilars are tuned to eliminate hub vibrations induced 
by forces referred to the rotating system at frequencies (N-1) Q 
or (N+1) Q, respectively. These forces are transformed to forces 
with frequency Nfi resolved to non-rotating axes. Thus, two separ- 
ate horizontal bifilars are required to eliminate Nfi forces re- 
solved to the fixed system. Post-processing of response results 
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obtained from the system solution may be required to evaluate 

horizontal bifilar designs. The Sample Applications section 

illustrates post-processing of y and y responses. A formula is 

o s 

given (see equation (101)) expressing the bifilar mass (N-1) Q 
displacement in terms of y^ and y^ to permit comparisons to be 
made with measured bifilar mass accelerations. 

Vertical Linear Bifilar BF2 

Vertical linear bifilar BF2 represents a bifilar vibration ab- 
sorber designed to reduce N per rev vibrations perpendicular to 
the plane of rotation of a helicopter rotor. The bifilar absorber 
is assumed to consist of N arms equal in number to the rotor 
blades. Each arm holds a pendular mass free to flap up or down. 
The arms are assumed to be equally disposed and each holds an 
identical absorber mass. The bifilar absorber rotates at the 
helicopter rotor speed Q. Figure 9 is a schematic one arm of the 
bifilar. The derivation of the equations of motion employed in 
the BF2 model is similar to the derivation described in the 
previous section for bifilar BFl. Each bifilar mass flapping 
displacement is expressed as the sum of a bifilar flapping collec- 
tive and cyclic displacements. 

^n " ^o^^ + PgSin i|i^) (2/N) (77) 

Substitution of equation (77) into the equations of motion (see 
previous section for a similar discussion) yields a nine co- 
ordinate constant coefficient equation set. The nine coordinates 

consist of bifilar mass collective (p ) sine cyclic (p ), and 

o s 

cosine cyclic (P^) amplitudes, and six hub displacements shown in 
Figure 9 . Restrictions on the configurations that can be analyzed 
are bifilars which have three or more equally disposed and identi- 
cal masses. Solution modes that may be used with BF2 are FRl, 
EG2 , and THl. Initial values of bifilar mass collective, and 
cyclic amplitudes may be specified for transient analysis. 
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N MASSES. ONLY ONE MASS SHOWN 
/8 =/3q/N + (i9c COStJ^+j8s S|N«|f)(2/N) 


Figure 9. Vertical Linear Bifilar BF2. 
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Vertical bifilars are tuned to eliminate hub vibrations induced by 
forces referred to the rotating system at a frequency equal to Nfi . 
These forces transform to Nfi forces resolved to non- rotating axes. 
Post-processing of response results obtained from the system 
solution may be required to evaluate vertical bifilar designs. 
The approach is analogous to that described for the horizontal 
bifilar (see previous section). 

Vibration Isolator ISl 

Vibration isolator ISl is a nodal isolator which may be placed 
between an excitation source, like the rotor/transmission system, 
and an airframe, to minimize vibrations on the airframe. The 
concept of the nodal isolator is described in Reference 7. The 
isolator employs a pivoted bar weight to balance spring and 
inertia forces to create a node or point of zero vibratory motion 
at the airframe attachment. When damping is present, vibrations 
are minimized at the airframe attachment at the excitation fre- 
quency. In coupled rotor/airframe systems the excitation fre- 
quency is N per rev where N is the number of rotor blades. Figure 
10 illustrates the basic nodal isolator, whose coordinates are 
upper and lower mass displacements. The antiresonant bar dis- 
placement is a function of upper and lower mass displacements. 

The upper mass simulates the rotor/transmission system. The lower 
mass simulates the airframe. Design variables are spring rate, 
antiresonant bar weight, bar weight overhang, and distance between 
bar pivots. Variables are selected to achieve minimum vibrations 
at a point on the airframe at a desired excitation frequency 
applied by rotor/transmission A. Design variables are selected to 
protect the isolator from resonance at the rotor excitation 
frequency. Reference 7 describes the dynamic attributes of the 
system. 
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Nodal isolator I SI generalizes the isolator concept to introduce 
six translational and rotational displacements at each upper and 
lower isolator attachment. The antiresonant bar weight is 
allowed to translate in three directions. Figure 11 illustrates 
the physical isolator and a schematic corresponding the physical 
isolator. Figure 11 does not show a second antiresonant bar 
weight included in the mathematical model. The lower bar weight 
permits additional design flexibility and is shown in Figure 12. 
The isolator transmission attachment (node 1) in Figure 11 corres- 
ponds to A in Figure 10. The isolator airframe attachment (node 
2) corresponds to B in Figure 10. Six coordinates at each attach- 
ment added to three coordinates for each antiresonant bar weight 
sum to a total of eighteen coordinates for ISl. Connection node 
numbers and Euler angles defining the isolator bar orientation are 
associated with each attachment. These attributes enable the 
isolator bars to be connected and assembled with any other sub- 
structures of the Base Program. A subsequent section (see Sample 
Applications) describes a configuration assembled from four 
isolator bars of type ISl to create a "six degree-of- freedom" 
nodal isolator. Nodal isolator ISl can be specialized to the 
isolator shown in Figure 10 . 

Rotor Aeroelastic Substructure RE3 and External Program G400 

Rotor aeroelastic substructure RE3 is used to represent a multi- 
blade aeroelastic rotor in forward flight or hover. The rotor 
representation in this substructure is defined by a rotor hub 
excitation force and moment vector and a rotor impedance matrix. 
The force and moment vector and impedance matrix may be generated 
by the G400 aeroelastic analysis. The Base Program combines the 
aeroelastic rotor representation with the impedance elements of 
other dynamical components to form a coupled system. The de- 
finitions of the rotor hub force excitation vector and rotor 
impedance matrix are given in a preceding section (Assembly of 
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Figure 11. Vibration Isolator ISl 
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Rotor Impedance Matrix). The coordinates represented in the RE3 
substructure are the six hub displacements listed in equation 
(59). The RES rotor impedance substructure is limited to the 
steady state forced response solution mode, FRl . The frequency of 
force excitation specified in FRl must be the same as that speci- 
fied for G400 to derive the rotor impedance matrix. This fre- 
quency must be an integer multiple of rotor speed and blade number 
(u)^ = nNfi). Component RES may be flagged to calculate Higher 
Harmonic Control (HHC) angles to minimize vibrations, such as 
described in a later section ( Higher Harmonic Control for 
Vibration Reduction). Optimal control angles, and blade stresses 
resulting from hub displacements and HHC control inputs may be 
calculated. This is achieved by utilizing data from G400 defining 
the effects on hub loads and blade stresses of HHC angle pertur- 
bations and hub displacements. Table 5 summarizes G400 
capabilities . 

Component RES may also be used to represent the impedance of a 
substructure which is not a rotor, such as an empennage. Aero- 
dynamic dampings may be represented in the impedance matrix. 
Structural and inertial contributions to the impedance matrix may 
also be included. 


HIGHER HARMONIC CONTROL FOR VIBRATION REDUCTION 


Higher harmonic control (HHC) is a concept for helicopter vibra- 
tion reduction which employs high frequency blade pitch controls 
to reduce helicopter vibrations. Vibration reduction is accom- 
plished by modification of vibratory forces and moments acting on 
the rotor blades, such that the resultant hub forces and moments 
are reduced. This section describes the approach used in the 
Coupled Rotor/Air frame Vibration Analysis to calculate helicopter 
vibrations for specified HHC angles and to calculate optimal HHC 
control angles for minimiim vibrations. 
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TABLE 5 - CAPABILITIES OF G400 ROTOR AEROELASTIC ANALYSIS 


Rotor Type 

Conventional articulated 
. Composite bearingless 
. Teetering 

Blade 

Redundant torque tube flex beam 
. Pendulum absorbers 
. Structural twist 
. Aerodynamic twist 

Aerodynamics 

Static aerodynamic data 
Unsteady data 

Nonuni form rotor induced inflow 

Stress Calculation 
Mode deflection 
. External force integration 

Mathematical 

. Uncoupled normal modes 
. Time history solution 

. Hub force and moment excitation vector and 
rotor impedance matrix 
. Higher harmonic control matrices 
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Higher harmonic control may be achieved by superimposing non- 
rotating swashplate motions at the blade passage frequency upon 
the collective and cyclic control inputs at fi. The blade passage 
frequency is NQ for an N bladed rotor. The Nfi cyclic motions of 
the non-rotating swashplate transform to (N-l)Q and (N+l)^^ blade 
pitch motions referred to the rotating system. The NQ collective 
motion of the swashplate results in NQ blade pitch changes in the 
rotating system. The harmonic coefficients of the blade pitch 
angle referred to the rotating system may be expressed as 


0 


T 

c 


0 


T 

S 


^®N-1 ' 

®N ' ®N+1}^ 

(78) 

II 

© 

®N ' ®N+1} 

(79) 


The elements of the vectors in equations (78) and (79) are the 
coefficients of cosine and sine terms in the pitch angle ex- 
pression with frequencies (N-l)Q, NQ , and (N+1)Q. The control 
inputs induce rotor forces and moments at NQ in the non-rotating 
system. 


Higher harmonic control angles are perturbed in the G400 analysis 

Matrix h^^ relates the rotor hub 
reaction forces to HHC control perturbations, and enables HHC to 
be introduced in the Coupled Rotor/Airframe Analysis. Successive 
perturbing of HHC controls enables columns of the transfer matrix 
to be identified. The procedure is similar to the determination 
of the rotor impedance matrix leading to equation (60). The 
resulting expression for the harmonics of hub reaction force at 
frequency NQ is 


to obtaxn a transfer matrix 
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( 81 ) 


"h ^ 

IIH * 12H 

- . » 

21H I 22H_ 

Equation (80) is the same as equation (60) with an additional term 

involving the HHC angles , Matrix h^j has dimensions 12x6 . Each 

submatrix ^i 2H' ^21H' ^22H equation (81) has dimension 

6x3. Each of the control vectors 6 and 6 has dimensions 3x1 as 

c s 

shown by equations (78) and (79). 

Assembly of the rotor hub excitation forces induced by HHC with 
the properties of other substructures is achieved by utilizing 
substructure RE3 . The assembly procedure parallels the hub force 
and moment assembly described in a preceding section (see Assembly 
of Rotor Impedance Matrix). The resulting system equation ex- 
pressed in terms of dependent coordinates is 


= 
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1 

‘^f^D ^12D 

-u)fCD + 

1 

1 
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<^D>s 




^IID : «12D 
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®21D » ^2 2D 


0 


e 


(82) 


Terms in equation (81) are the same as the terms in equation (64) 
except for the HHC terms. Matrix in equation (82) is 
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( 83 ) 


The matrix in equation (83) is zero except for 6x3 element 
which derives from equations (80) and (81). The forms of matrices 
^21D^ ^22D equation (82) are similar to equation (83). 

The row dimensions of the zero rectangular matrix [0]^^^ is equal 
to the numher of coordinates for the i— th substructure in the 
system. The column dimension of this submatrix is 3. 


The transformation from dependent to independent coordinates for 
the system is analogous to the derivation in equations (67) to 
(73). The resulting system response is 



(84) 


X 


CO 


X 


so 



(85) 
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(89) 


Equation (85) is the same as equation (27) derived in a preceding 
section (Assembly of Rotor Impedance Matrix). Equation (85) 
yields the system response including a rotor, without HHC. 
Equation (84) yields the system response when HHC is specified. 
It is essentially the same as equation (1) in Reference 11. 
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Matrix T is the transfer matrix relating HHC inputs to vibratory 
response. The dimensions of the matrices introduced in equation 
(84) are: T = 2Nx6, A = 2Nx2N, H = 2Nx6, ®21' *^22 “ 

where N is the number of coordinates in the independent coordinate 
system. 

To derive the control angle vector for minimum vibrations , a 
scalar performance index J is defined. 



The control vector must be chosen to minimize J. The vector of 
response variables Y, whose cosine and sine components are and 
Y^, may be identified with the displacement vector the inter- 
face force vector or a vector composed of both displacements 

and interface forces. In this application, Y is selected as 

Y = (91) 

which identifies Y as the dependent coordinate vector. The 

diagonal weighting matrix W permits the weighting of some 

z 

responses more than others . Diagonal elements may be selected to 
be zero except at a point on the structure whose displacements are 
to be minimized. A typical requirement is the weighting of 
cockpit or tail responses for minimization. Weighting matrix 
has dimensions x where Nj^ is the size of the vector of 

dependent coordinates, 
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The 3x3 weighting matrix limits the amplitudes of the rotor 

controls. Equation (90) is analogous to equation (2) in Reference 
11 . The control state for a minimum value of J is obtained by 
first substituting equation (84) into equations (67) and (68). 
This yields 



Equations (91) and (92) are substituted into equation (90) to 
express the performance index in terms of control variables. The 
resulting expression is differentiated with respect to the control 
variables and equated to zero to obtain the control state for a 
minimum value of J. 



^ W p J 0 
2 I 

1 

0 ; 


T 


+ 



0 

T^ 

r p^w p 1 0 “1 

zP , 

L 0 

1 

CD 

> 

1 

L 0 1 P W^p J 



(93) 


Equation (93) is analogous to equation (3) in Reference 11, except 
that it is expressed in terms of cosine and sine components of 
the independent coordinate vector, . An equivalent version of 
equation (93) is implemented in the program code to minimize the 
computer storage requirement. Transformations are applied to the 
T matrix in the code as opposed to the matrix such as shown in 
equation (93). 
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Substituting equation (93) into equation (92) yields the minimum 
vibration vector. The optimal control solution and corresponding 
vibrations are calculated in the Base Program. File data are read 
to acquire the transfer matrix h^ tramsmitted from G400. A 
non-linear dependence of h^ on the control state is approximated 
by deriving h^ for control perturbations centered at prior solu- 
tions obtained from the Base Program. The HHC solution is refined 
by iterating between the Base Program and G400 . 


Blade stresses and push rod loads are monitored for the effects of 
rotor hub displacement and higher harmonic control inputs. This 
is achieved by expanding the harmonic coefficient of a stress in 
terms of hub displacements and HHC control angles, such as was 
done for the harmonics of hub reaction force in equation (80). 


n 


= (S ) 

' n^o 




i c H 

f X „ 

\ s H 



(94) 


Stress coefficient S is the n-th harmonic of stress at a point on 

n 

the rotor blade where n is an integer multiple of rotor speed, Q. 

Stress coefficient S^ represents either the flatwise, edgewise, 

torsion stress, or push rod load. Stress harmonic coefficient S^ 

is a scalar. The coefficient (S ) is the n-th harmonic coeffic- 

n o 

ient of stress without hub displacement nor HHC controls. Row 

matrix ^ 1x12 matrix of coefficients introducing the 

effects of hub displacement. Row matrix is a 1x6 matrix of 

coefficients introducing the effects of HHC controls. Scalar 

(S ) and matrices [S and [S _] are calculated by G400 and are 
'no nH n0 

transmitted to the Base Program. The Base Program calculates the 
harmonics of stress or push rod load from equation (94). These 
harmonics are used to calculate the variation of stress or push 
rod load with time. Peak to peak values are determined in the 
Base Program. 
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SAMPLE APPLICATIONS 


A discussion is given below of sample applications performed with 
the Base Program to demonstrate the capabilities of the analysis 
package. The set of applications exercises each component 
(substructure or solution mode). The problems were selected for 
physical significance and ability to test program capabilities and 
component combinations . Table 6 summarizes the attributes of the 
Base Program sample problems. A brief description is keyed to a 
list of components employed for the case. 

Base Program Applications 

Case 1 was used to verify the functioning of the horizontal linear 
bifilar BFl and components GFl, MSI and FRl . Attributes of the 
UH-60 system were used. The configuration consists of a bifilar 
with four pendulum masses, coupled to a modal fuselage with five 
modes. Shake test data vjere the source of modal properties- A 4Q 
inplane force is applied to the hub to simulate the effects of a 
rotor. The components of the applied force are resolved to the 
fixed system and are equal to 




^xl = 

1700 

cos 

4f2t, 

N (lb) 


(95) 




1100 

sin 


N (lb) 


(96) 

The above 

values of 

the 

force coefficients are 

input 

to GFl . 

Figure 

2 

shows the 

positive 

directions of the 

forces . 

These 

forces 

are 

equivalent 

to 








^X2 = 

1400 

cos 

3fit 

+ 300 cos 5fit, N 

(lb) 

(97) 



II 

1400 

sin 

3Qt 

- 300 sin 5Qt, N 

(lb) 

(98) 
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excitation force. 

Elastic beam model comprising V V v 7840 
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TABLE 6 . - Continued 















referred to the rotating system. The fixed system 4Q forces 
result from rotating system forces with frequencies and 5Q 

referred to the rotating system. The rotating system forces 
applied by the rotor were found from measured bifilar arm ac- 
celerometer readings. Hub displacements were assumed to be 
negligible as a result of the bifilar having minimized the hub 
motions. The steady state forced response of the system was 
obtained from solution mode FRl. The coupled system solution has 
cosine and sine components (see equation (24)) and the bifilar arm 
displacements are therefore 

(Yo)cs 

S S o o 

Y = (Y ) cos 4fit + (y ) sin 4Qt, rad (100) 

c c c c s 

Substitution of equations (99) and (100) in equation (74) yields 

A/an = i + [(Vs>s + rad (101) 

4 


For the bifilar mass 3Q amplitude additional cases were run to 
find the bifilar mass amplitude variation with airspeed, shown in 
Figure 13 . Results are compared with test data obtained from a 
bifilar mass mounted accelerometer. The agreement of theoretical 
results and test data shows the correctness of the substructure 
coupling procedure, the equations of motion for the components 
employed in the configuration, and the system properties deter- 
mined from test data. The analytical model is not suitable for 
bifilar or fixed system design when rotor hub loads data are not 
available or at bifilar off-design conditions with large hub 
displacements when rotor impedance effects are important. A 
suitable general model for design studies would employ an aero- 
elastic rotor model, such as RE3 . 
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Figure 13. Bifilar Mass Amplitude Variation with Airspeed 


Case 2 is a verification of a vertical linear bifilar BF2. A 
vertical rotor induced force is applied to the hub. In other 
respects, the case is similar to Case 1. 

Case 3 is a verification of beam model BMl, constraint component 
CNl , and eigensolution mode EGl . The case yields the mode shapes 
and frequencies for a five element uniform beam. The beam is 
constrained from motion in all directions except vertical trans- 
lation. Figure 14 shows mode shapes and frequencies for a ten 
element beam model similar to Case 3. Results are in close 
agreement with classical results. The ten element beam case 
required approximately ISK^^ decimal working storage words (IK^q = 
1024 words ) , which together with the storage required for instruc- 
tions and common blocks (IVK^^ decimal words) shows that this 
problem has approached the 32K^^ decimal word limit for inter- 
active CDC system usage (see Reference 17). The beam model BMl 
could have been used with any solution mode. In particular, a 
forced response study with FRl and rotor system RE3 could have 
been performed to study the vibration behavior of a helicopter 
preliminary design - 

Case 4 is a verification of the fixed absorber FAl, with parameter 
variation mode PVl . A mass is attached to the base of the ab- 
sorber by means of an MSI component to simulate a two mass system 
joined by a spring and damper. An oscillating force is applied to 
the base of the absorber and a steady state forced response 
solution is obtained. The forcing frequency is specified to range 
between one and ten Hz. This parameter variation is achieved by 
means of input to PVl. Vibratory response results agree with 
analytical derivations. Analytical results were derived manually 
by application of the Laplace transform method to the assumed 
system. 
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Case 5 is similar to Case 4 except that a time-history solution 
mode THl substitutes for steady-state forced response mode FRl and 
parameter variation mode PVl, The case is also used to verify the 
general (GEN) control mode line printer output suppression capa- 
bility. The case transmits data to a plot file and the Plot 
Package is used to plot results. Figure 15 illustrates the 
variation with time of the responses of the two masses. The 
steady-state values of vibratory response agree with the values 
obtained from Case 4, and with analytical derivation. A rigid 
body component of displacement of the two mass system shows a 
linear increase with time and agrees with analytical derivation. 

Case 6 is the same as Case 5 except that it is used to check the 
restart option in the time-history solution mode. The case 
employs time-history restart data transmitted to a file from a 
prior time history run. 

Case 7 consists of one nodal isolator of type I SI coupled to upper 
and lower masses, connected by a spring such as shown in Figure 
10. The isolator is free to move vertically. The frequency 
response of the system to a force applied to the isolator upper 
node is shown in Figure 16. The isolator capability to create a 
point of zero vibration at the lower attachment node is evident at 
the antiresonant frequency (27 Hz). 

Cases 8 to 10 serve as a systematic synthesis and verification of 
features employed in a hypothetical six degree-of- freedom 
isolator model. Case 11. 

Case 8 is similar to Case 7 in the character except that the 
isolator is free to move vertically and laterally. Two modal 
structures (MSI) were connected to the isolator at nodes one and 
two. A force of type GFl was applied to node 1 of the isolator. 
The configuration is shown in Figure 17. 
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Figure 16. Frequency Response of a Simple Nodal Isolator. 
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Figure 17- Configurations Used to Test I SI Component. 






Case 9 is similar to case 8 except that the GFl force was applied 
to node 3 of an MSI substructure (see Figure 17) which connects 
the force application point and the substructure attachment to the 
upper isolator node (node 2). The upper MSI substructure simu- 
lates a transmission placed between a rotor and an isolator. The 
lower MSI substructure simulates a fuselage. 

Case 10 generalizes Case 9 to a configuration employing six modes 
for transmission representation, six modes for fuselage representa- 
tion, and an isolator of type ISl. All the modes are rigid-body 
modes. A steady state forced response solution is obtained. 

Case 11 consists of a six degree-of- freedom isolator constructed 
from four isolator bars of type ISl connected at right angles to 
each other. Six rigid-body transmission modes, and six rigid-body 
fuselage modes are simulated. Six hub forces and moments are 
applied to the upper transmission to simulate the loading applied 
by a rotor. Isolator characteristics were derived from a half 
scale UH-60 isolator. An illustration of the model is shown in 
Figure 18. A total of nine connection nodes is used. This case 
introduces substructure Euler angles to check the correct forma- 
tion of angular transformations from local to global axes. Bar 
Euler angles of 0, 90, 180 and 270 degrees are introduced. Figure 
19 shows the frequency response of the isolator at a fuselage 
attachment. The variations with forcing frequency of three 
components of translation are shown at this point. Figure 20 is a 
three dimensional view of vertical displacements at the four 
fuselage attachments (see Reference 18 for plot types that can be 
requested) . 

Case 12 verifies the ability of the program to employ rotor 
aeroelastic model RE2 (E927 rotor) . A five mode fixed system was 
coupled to the rotor. The rotor employs twelve coordinates 
consisting of flap collective, lead collective, two flap cyclic 
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Figure 18. Six Degree-of-Freedom Isolator, 



AMPLITUC 



FREQUENCY (HZ) 


Figure 19. Frequency Response of Six Degree-of -Freedom 
Isolator at a Fuselage Attachment, 



AMPLITl 



Figure 20. Six Degree— of— Freedom Isolator Vertical 
Responses at Fuselage Attachment Points. 



modes, two lead cyclic modes, two bending collective modes, and 
four bending cyclic modes. An eigensolution of the system was ob- 
tained using the eigensolution mode with damping, EG2 . Results 
agreed with results obtained from the unmodified E927 analysis, 
which included an elastic fuselage. 

Case 13 verifies the ability of the program to employ rotor 
aeroelastic analysis RES and to serve as a tool to study effici- 
ently the effects of variations in components internal to the Base 
Program. A rotor hub load excitation vector and a hub impedance 
matrix were transmitted by G400. A hypothetical full scale four 
blade rotor and four fixed system modes were specified. One flap, 
one lead, three uncoupled flatwise, two uncoupled edgewise, and 
one uncoupled torsion blade modes were used. Airspeed of 120 
knots was specified, and uniform inflow was used. Figures 21 and 
22 show the effects of changes in the frequency of the third 
fuselage mode on hub loads and accelerations, respectively (magni- 
tudes of loads and accelerations should not be interpreted as 
being representative of an actual helicopter). Selected elements 
of the rotor impedance matrix were derived analytically and agreed 
with the G400 results. Verification of the ability of program 
package to derive and use rotor hub loads and impedances for 
rotors employing variable inflow is covered in the subsequent 
text. Case 13 verifies the END CASE command (see Reference 17). 
The first set of data in Case 13 employs a flag in RE3 to display 
the output solution vector in feet and the rotor hub interface 
forces (see equation (60)). The data following the END CASE 
command request the output to be displayed in units of accelera- 
tion of gravity (see equations (28) and (29)). 

EXTERNAL PROGRAMS IN THE PROGRAM PACKAGE 


Descriptions are given here of the functions and assumptions of 
the External Programs available in the Coupled Rotor/Airframe 
Vibration Analysis package. 
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Figure 22. Effects on Rotor Hub Accelerations of Variations 
in Third Fixed System Mode Frequency. 
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Rotor Aeroelastic Analysis 'E927 


The E927 Rotor Aeroelastic Stability Analysis (see Reference 6) is 
a linear analysis yielding dynamic, and aerodynamic mass, damping, 
and stiffness matrices of the rotor system in hover. These 
matrices are associated with perturbations of the coordinates of 
the rotor system. Coordinates comprise blade generalized dis- 
placements, both physical (e.g. flap angle) and modal. Five hub 
displacements are represented; lateral, longitudinal, and vertical 
translations, and hub pitch and roll. Linearization occurs about 
a rotor trim state. The trim vector is obtained from G400 or any 
other suitable analysis. Mass, damping, and stiffness matrices 
are transmitted to a file and read by the Base Program. These are 
assembled with matrices for other components selected. The Base 
Program eigensolution mode is used to evaluate the stability of 
the rotor system coupled to other components . Steady-state forced 
response may be calculated if a rotor hub force and moment excita- 
tion vector is input to the Base Program. Hub force and moment 
excitation data are obtained from the G400 analysis or any other 
suitable analysis. 


Rotor Aeroelastic Analysis G400 

The G400 Rotor Aeroelastic Analysis yields the time-history 
response of conventional articulated, teetering, and composite 
bearingless rotors. Options permit the use of variable rotor 
induced inflow, static blade aerodynamic data, and unsteady 
aerodynamic data. Blade mounted pendulxim absorbers are modeled. 
Forward flight and hover are represented. The program calculates 
the hub excitation vector for a fixed hub, and the rotor impedance 
matrix relating six hub forces and moments to harmonic perturba- 
tions of the hub. These data are transmitted to a file which is 
read by the Base Program. The rotor hub excitation vector and the 
hub impedance are assembled with mass, damping, and stiffness 
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elements of internal components in the Base Program. The steady- 
state forced response mode in the Base Program is invoked to 
obtain coupled system responses . 

Rotor aeroelastic analysis G400 was created to overcome inadequa- 
cies of applications of conventional analyses to composite bearing- 
less rotors. The analysis is based on a noonal modes model of a 
non-linear rotor blade. References 8 and 9 describe -the basis of 
the method. A coordinate transformation is applied to second 

derivatives of blade bending deflections as opposed to bending 
deflections. This transformation enables uncoupled normal modes 
to be used to achieve non-linear twist representations which 
satisfy the blade root and tip boundary conditions. Uncoupled 
modes derive from eigensolutions which employ uncoupled flatwise 
bending, edgewise bending, and torsion equations. Up to four 
flatwise, three edgewise, and three torsional elastic modes may be 
used. Blade articulation reduces the number of allowed modes. A 
Galerkin method yields modal differential equations which are 
integrated to obtain the time history response of the blade. 
Blade stress calculation options are mode deflection and external 
force integration methods (see Reference 10). Single blade loads 
are analyzed harmonically to obtain harmonics of rotor loads for a 
fixed hub. The hub is perturbed by specified harmonic motions and 
the resulting hub load perturbations are analyzed harmonically to 
obtain the rotor impedance matrix (see Assembly of Rotor Impedance 
Matrix). Control perturbations yield matrices required with the 
HHC solution mode (see Higher Harmonic Control for Vibration 
Reduction). Table 5 lists the capabilities of the G400 analysis. 
All rotor types are based on a single blade analysis, except for 
teetering rotors which are treated as two bladed rotors. Variable 
rotor induced inflow is introduced by executing sequentially G400 
and the F389 Prescribed Wake Rotor Inflow Analysis. Refinement 
with variable inflow precedes calculation of the hiib excitation 
vector and impedance matrix. The same interfacing of the Base 
Program and G400 occurs as in the uniform inflow case. Fuselage 
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induced inflow at the rotor disc is introduced by transferring 
harmonics of these data from WABAT to inflow program F389. 

Prescribed Wake Rotor Inflow and Flow Field Analysis F389 

Prescribed wake rotor inflow and flow field analysis F389 cal- 
culates the aerodynamic flow field induced by a rotor. The 
analysis may be used to calculate rotor induced velocities in the 
plane of the rotor or at user selected points - e.g., at the 
empennage . 

The F389 analysis is based on a lifting line rotor blade re- 
presentation. A vortex wake of prescribed shape is assumed to 
trail from the rotor blades. The assumption built into the 
program is that the wake form is a skewed helix (classical wake) 
as shown in Figure 23 . The skew angle and distance between 
successive helixes is calculated by assuming the wake to be 
convected from the rotor by a uniform self-induced rotor inflow 
derived from momentum balance considerations, and from contribu- 
tions from components of vehicle flight speed. The distance from 
the blade at which the trailing vortices roll up into a single tip 
vortex may be specified. The wake is terminated at a finite 
number of wake revolutions. The method of finding unknown circu- 
lations on blades and wake is described in Reference 12. The 
state of motion of the blades is assumed to be known. This state 
is derived from any suitable analysis external to F389- 

The procedure consists of using the Kutta-Joukowski condition to 
derive a formula expressing blade circulation as the product of 
blade section lift curve slope, and blade section effective angle 
of attack. The effective angle is the geometric angle of attack 
corrected for inflow induced by rotor blades and wake. Circula- 
tion continuity (Helmholtz law) is used to express the circulation 
values in the trailing wake in terms of the circulations on the 
blades. A steady-state of flight is assumed, and as a result, the 
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23 . 


Computer Wake Representation-Undistorted Wake 


Figure 





azimuthal variation of circulation is periodic. Flow velocities 
induced by blades and trailing vortex wake are found from the 
Biot-Savart law and contain circulation variables as unknowns. 
The induced velocity expressions are substituted in the formula 
for blade circulation to obtain a finite set of linear simultan- 
eous equations which are solved for the unknown circulations . 

After the circulations are known, induced velocities can be 
calculated anywhere in the flow field by application of the Biot 
Savart law. 

Section airfoil lift coefficients are approximated by straight 
line curve fits to airfoil data, approximating unstalled and 
stalled regions of lift curves (see Figure 24). 

The F389 program can be made to calculate the flow field consis- 
tent with the motions of an elastic blade. This is achieved by 
calculating the component of airflow normal to the median surface 
of the rotor blade induced by the motion of the elastic blade, 
which excludes the rotor induced flow contribution to the airflow. 
This airflow component is 

w, = U„0 + U - ROTOR INDUCED FLOW (102) 
1 T p 

The calculation of w^ is performed in G400. The angle, 0, is the 
blade pitch angle. The variable, w^, is automatically transmitted 
to F389 by G400. Program F389 adds the rotor induced inflow to w^ 
to calculate the effective angle of attack of a blade section. 
Program F389 can represent an elastic rotor since no assumptions 
are made restricting and to non-elastic blades. It should 
be noted that another mode of F389 operation is to employ the 
flapping harmonics to define the blade motion input. This mode of 
calculation is less accurate than the mode employing w^ because it 
omits elastic displacement contributions to airflow velocity 
components . 
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START 



STOP 


N (Execution 
Flag) 


Execution Sequence 


-2 

-1 

0 

1 

2 


F389 only 
G400 only 
G400 only 

G400/F389/G400 (one cycle) 
G400/F389/G400/F389/G400 (2 


Figure 25. Calculation Sequences For The Coupled 

G400/F389 Programs 


cycles ) 
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Determination of rotor induced flow requires a coupling of F389 
and a rotor aeroelastic analysis. In the vibration analysis 
package, this is achieved by transmitting data automatically 
between F389 and G400. Data transmission and the invoking of F389 
and G400 is controlled by a Fortran driver program, VIBIT. Driver 
program VIBIT causes G400 to be run first. Uniform rotor induced 
inflow is used. Program G400 transmits w^ to F389 which calcu- 
lates the rotor induced flow. The induced flow is transmitted by 
F389 to G400 which recalculates the blade aeroelastic response. 
The execution sequence G400/F389/G400 is referred to as one cycle 
of operation. An arbitrary number of cycles may be performed 
until the user is satisfied that consistency of blade circulations 
(lifts), blade motions, and aerodynamic velocities is achieved 
among the programs. Figure 25 illustrates the VIBIT/G400/F389 
calculation sequences . After a converged solution is achieved 
calculations may be performed in the same run to obtain rotor 
output data - e.g., hub excitations and impedance. 

The F389 program may be used to calculate rotor induced velocities 
at the empennage to provide data for the calculation of rotor wake 
induced empennage excitations by the RIEVA program (see Figure 1). 
The empennage flow field is assumed to have no effect on the rotor 
and no iterative calculations involving F389 and RIEVA are per- 
formed. 

The F389 program accepts input defining velocities induced by a 
fuselage or gust. The velocity field is defined as harmonics of 
inflow at the rotor. Fuselage flow field data may derive from 
WABAT (see Figure 1). No change in fuselage velocities is assumed 
to occur as a result of the rotor wake influence on the fuselage. 

Wing and Body Aerodynamic Technique (WABAT) 

The Wing and Body Aerodynamic Technique (WABAT) is used to calcu- 
late a fuselage flow field. The flow velocities may be obtained 
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on the surface of the fuselage or at arbitrary (off-body) points. 
The velocity components in the plane of the rotor can be expressed 
as a harmonic series of blade azimuth angle whose coefficients 
depend on blade radius. The harmonic coefficients are transmitted 
to F389 to define the effects of fuselage induced flow at the 
rotor . 

To obtain the solution to fuselage flow, the surface of the 
fuselage is divided into planar quadrilateral panels on which 
aerodynamic sources of constant strength are placed. The flow is 
assumed to be incompressible - i.e., the airspeed of the heli- 
copter is assumed to be low. The sources satisfy the Laplace 
equation, which governs incompressible, inviscid (potential) flow. 
Strengths of the panel sources are determined from the condition 
that the flow normal to the fuselage is zero or that it is some 
specified value like that induced by an engine intake or exhaust. 
The number of solution points is specified to be equal to the 
number of sources of unknown strength. A set of linear simultan- 
eous equations for the panel source strengths is obtained and 
solved. When the source strengths are found, the velocity field 
induced by the fuselage at any field point may be calculated by 
summing the contributions of the panels. Appendix B of Reference 
13 gives further details of the method. 

Lifting wings are represented by distributing horseshoe vortices 
(see Figure 26) in the mean surfaces of the wings and in the wing 
wake (trailing vortices). The solution technique for wing and 
body configurations consists of solving for the distribution of 
lift on the wing without the body. This lifting solution is 
achieved by means of a modified Multhropp analysis (see Appendix B 
Reference 13 ) . The Multhropp lift distribution assigns circu- 
lation weighting factors among the horseshoe vortices. In other 
words, the relative strengths of the circulations among the 
horseshoe vortices distributed on a chord at each spanwise station 
are known. The magnitudes of the circulations are not known. As 
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a result of assignment of circulation weighting factors and a 
summing of the circulations on each chord, the circulation of each 
vortex can be expressed as some fraction of the total circulation 
at a spanwise station. The number of unknowns associated with the 
lifting wing at each spanwise station is therefore reduced to one, 
which is the total circulation at the station. The total number 
of lifting wing unknowns is N, where N is the nximber of spanwise 
stations on the wing. The determination of the circulation 
strength magnitudes is preceded by placing source panels on the 
wing to represent the thickness effects of the wing on the flow. 
The unknown circulations and source strengths are then obtained by 
summing the induced velocity effects of body and wing source 
panels and wing horseshoe vortices, and equating these velocities 
to specified velocities on the wing and body surfaces. When a 
wing is present, solution points are defined at each panel point 
on wing and body. In addition there are N solution points placed 
at the 3/4 chord of the aft-most horseshoe vortex at each spanwise 
station. The N additional points are introduced to increase the 
number of boundary conditions to accommodate -he N additional 
unknowns (circulations) of the lifting wing. As a result, the 
number of unknowns and solution points are equal and the equations 
are solved for source and horseshoe vortex strengths . A wing and 
body solution satisfying the boundary conditions of the fuselage 
and wing is consistently obtained. 

The synthesis of the fuselage geometry is facilitated by an 
automated technique of cross-section development from input data. 
Input data define a limited number of points on the cross section. 
This is performed by the Automated Paneling Technique (APT/Plot) 
program. Appendix A, Reference 13 describes the basis of the 
method. Incorporated in this program is a Tektronix plot capa- 
bility which permits the user to verify the correctness of the 
panel input and transmits the panel data by means of a file to the 
WABAT potential flow program. The user may develop the fuselage 
geometry in subsections . Each section consists of a number of 
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panels . The Tektronix plot package assembles the sections into a 
complete fuselage and displays the resulting fuselage. The 
division of the fuselage into sections facilitates the development 
of the complete fuselage geometry. Parts of the fuselage - e.g., 
engine nacelles - requiring complex panel geometries may be 
developed and verified separately. The plot package is able to 
provide views of the fuselage from any angle and distance to 
permit the user to study and understand the data. Sections of the 
fuselage may be separated from each other to understand the data 
and character of the sections. Figures 27 to 29 illustrate views 
of the fuselage derived from the APT/Plot program. Reference 14 
provides user instructions for the APT/Plot program. 

Panel geometry data are transmitted by means of a file to the 
WABAT potential flow program (Figure 1) to calculate the wing and 
fuselage induced flow. The matrix of aerodynamic influence 
coefficients of the wing and body singularities (sources and 
vortices) may be saved on a file for subsequent applications. 
Employment of a stored influence coefficient matrix reduces the 
run time for subsequent cases using the same wing/ fuselage geo- 
metry when attitude ( angle-of-attack, sideslip) is varied. User 
instructions are provided in Reference 14. A secondary mode of 
the APT/Plot package may be invoked to obtain a fuselage surface 
flow display after the WABAT program has been run. The arrows in 
Figures 27 to 29 illustrate the flow directions of perturbational 
flows on fuselage panels. 

Rotor Wake Induced Empennage Vibrations (RIEVA) 

The Rotor Induced Empennage Vibration Airloads (RIEVA) program 
calculates the harmonics of empennage lift, rolling moment and 
pitching moment induced by a rotor wake at the blade passage 
frequency (NQ). These harmonic coefficients are input to the 
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Figure 27. APT/Plot Program Fuselage View from Below- 



Figure 28. APT/Plot Program Fuselage View from Above. 
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Coupled Rotor/Airframe Vibration Analysis Base Program to calcu- 
late the vibratory response of a helicopter to rotor wake excita- 
tions of the empennage. The steady state forced response mode in 
the Base Program is used to obtain a solution for the coupled 
system . 

The empennage is represented as a three dimensional lifting 
surface of arbitrary shape. The empennage may include horizontal 
and vertical stabilizers. The flow is assumed to be incompres- 
sible and unsteady. The wing and its wake are represented by a 
doublet lattice on the mean surface of the wing and a vortex 
lattice on the wake. Figure 30 shows the lattice geometry. The 
perturbation velocities induced by the doublet and vortex lattices 
are calculated by the Biot-Savart law. At each instant of time 
the strengths of doublets and vortices must satisfy the boundary 
conditions of no flow normal to the wing surface, the Kutta 
condition of smooth flow from the trailing edge, and continuity 
relationships among the wake vortices and at the trailing edge 
interface between doublets and vortices. The distribution of 
circulations among the trailing vortices is known because vortices 
of known strengths are convected from the trailing edge at the 
speed of the flow. At the trailing edge there is an equivalence 
between the lifting surface doublet strengths and the strengths of 
the shed vortices. As a result of these relationships the number 
of unknowns reduces to the number of constant doublet panels on 
the empennage- These unknowns are assumed to satisfy as many 
simultaneous equations for normal velocities as there are boundary 
condition points specified on the empennage. The doublet 
strengths are obtained by simultaneously satisfying these equa- 
tions at each time instant. The wake is assumed to begin at the 
calculation starting time. When doublet strengths are known 
empennage surface pressures may be calculated at each instant of 
time. The passage over the empennage of rotor vortices causes the 
flow over the empennage to be unsteady. When the empennage 
aerodynamic response is determined to be periodic at the blade 
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Figure 30. Empennage and Wake Lattice Geometry in RIEVA Analysis 


passage frequency (Nf^) the time -stepping calculation is termin- 
ated. Surface pressures are integrated to obtain empennage lift, 
pitching moment and rolling moments. These are harmonically 
analyzed at frequency NQ . Vibratory response of a helicopter to 
rotor wake induced empennage excitations is obtained by input of 
these harmonics to the Base Program. Empennage loads are assumed 
to be applied to this fuselage junction. 

Close passage of rotor induced vortices to the empennage is 
approximated by assuming surface pressures to be induced by Ran 
kine viscous vortices (see Reference 15). These consist of 
vortices with a solid core and a potential vortex region outside 
the core. The vortex is assumed to be more than one vortex core 
radius away from the lifting surface. If the vortex is closer 
than one core radius to the lifting surface, it is moved to a new 
position such as described in Reference 15. 

The rotor induced wake is divided into a near wake region and a 
far wake region such as shown in Figure 31. The near wake region 
is contained by two vertical streamwise planes parallel to the 
flight direction, which also contain the empennage. The far wake 
is outside this region. The RI.EVA program provides an internal 
calculation of the near wake from data defining the flight speed, 
tip path plane angle, position of the empennage, and rotor momen- 
tum induced flow. The far wake effects are represented as rotor 
induced flow velocity harmonics . These harmonics are obtained 
from F389 and are input into RIEVA. The user may elect not to 
input far wake harmonics if these effects are small. (Omission of 
far wake data simplifies RIEVA program usage, and has been found 
to be justified on the basis of applications made to this date.) 

The following approximations are made in the present version of 
the RIEVA analysis. It is assumed first that the air pressures 
acting on the empennage are not influenced by the motion of the 
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Figure 31. An Example of Splitting of Rotor Wake into Two Parts 
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empennage. A second approximation is that the empennage dis- 
placement consists of plunge, pitch, and roll - i.e., the empen- 
nage is assumed to displace as a rigid body. The corresponding 
generalized forces acting on the empennage are lift, pitching 
moment, and rolling moment. These are assumed to act on the 
elastic structure at the empennage fuselage junction. When a 
modal representation of the fuselage is used it is necessary to 
define the components of modal displacement at this junction - 
e.g., plunge, pitch, and rolling elements of the modal matrix - to 
calculate the generalized forces induced by the empennage. 
Refinement of the calculation of generalized forces can be accomp- 
lished by multiplying empennage elastic mode shapes by the dis- 
played values of pressures on the lifting surfaces to obtain 
generalized forces. The code does not provide this calculation. 

The assumption that air pressures on the empennage are not in- 
fluenced by motion of the empennage is equivalent to neglecting 
aerodynamic damping in plunge, pitch, and roll, and aerodynamic 
stiffness in pitch. These effects may be introduced as follows. 
Aerodynamically induced impedance elements can be calculated by 
employing an unsteady aerodynamic analysis. A simple approach 
would be to use Theodorsen unsteady aerodynamics (see Reference 
16) and a strip theory assumption. Introduction of the empennage 
impedance matrix may be achieved by connecting a substructure of 
type RE3 to the point of application of empennage forces (see the 
subsection Rotor Aeroelastic Substructure RE3 and External Program 
G400 for use of RE3 as a general impedance substructure). 

There are several alternative ways to representing the structural 
and inertial properties of the empennage. One method is to 
provide a modal representation. Fuselage modes would be deter- 
mined with the empennage assiimed to be present. An alternative 
approach is to calculate fuselage modes without the empennage. 
Empennage properties would be represented as contributions to the 
impedance matrix of substructure RE3 , which are connected to the 
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fuselage at: the empeimage/fuselage junction. Impedance elements 
would also include aerodynamic contributions . 

A third approximation in the analysis is that it does not include 
the contribution to the streamwise velocity component induced by 
vorticity shed from the rotor. This increment may be important 
when a rotor vortex passes near the vertical stabilizer. Com- 
ponents of rotor induced velocities normal to the lifting surfaces 
and perturbational streamwise components induced by the lifting 
surfaces are represented and are essential for obtaining the RIEVA 
solution. 


Base Program Plot Package 

The Base Program stores solution information on file for pro- 
cessing by the Base Program Plot Package. Line printer output 
produced by the Base Program may be suppressed optionally. The 
Plot Package is designed to recognize and display solution modes 
(forced response, time history, eigensolution) and coordinate 
labels without user intervention beyond identification of co- 
ordinate numbers selected for plotting. Comparative and per- 
spective plots are options (see Reference 18). 

EXTERNAL PROGRAM APPLICATIONS 

Table 7 lists the check cases used to verify the capabilities of 
the External Programs. There are six external programs comprising 
the linked G400/F389, Tail Option (RIEVA), modified E927, Base 
Program plot package, WABAT geometry program/plot package (Auto- 
matic Paneling Technigue/Plot or APT/Plot), and WABAT potential 
flow program. The check cases were selected to verify key capa- 
bilities, the abilities of the programs to transmit data to the 
Base Program and to each other, and to acquire data from the Base 
Program, such as plot data. 
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TABLE 7 - INPUT CHECK CASES FOR EXTERNAL PROGRAMS 
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Coupled Rotor Aeroelastic Analysis and 
Variable Inflow Analysis G400/F389 


Case 1 verifies the ability of the G400 module to calculate an 
articulated rotor response. A full scale rotor of 22 foot radius 
was specified. Blade properties are typical of actual rotors. 
Uniform inflow was specified. This case uses specified input 
control angles - i.e., trim was not requested. Blade coordinates 
are flap, lead, three elastic flatwise, two edgewise, and one 
torsion uncoupled modes. 

Case 2 is similar to Case 1, except the blade pendulum vibration 
absorbers are specified. 

Case 3 is an application of G400 which is used to check the 
formation of the hub load excitation vector, rotor impedance 
matrix, matrix h^^ containing the effects of input HHC angles on 
hub loads and blade stress matrices (see previous sections for 
descriptions of the above terms ) . The data are transmitted to a 
file- Base Program check case 13 (Table 6) uses the hub load 
excitation vector and rotor impedance data to calculate system 
response without HHC. Numerical rotor trim is requested, as 
opposed to analytic (Bailey derivative) trim (see Reference 9). 
Analytic trim was invoked for the HHC ARES model applications 
described subsequently. In other respects, the case is similar to 
Case 1 . 

Case 4 is similar to Case 3 except that a teetering rotor is 
specified. 

Case 5 exercises the F389 program in a stand-alone mode. Inflow 
harmonics data are transmitted to a file which may be read by the 
Tail Option (RIEVA) module to represent the effects of the rotor 
induced far wake at the tail. 
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Case 6 exercises the F389 program in a stand-alone mode. Inflow 
harmonics data are transmitted to a file by F389 which may be read 
by G400, employed in a stand-alone mode. 

Case 7 is a coupled G400/F389 run. Two cycles of calculation 
between G400 and F389 are performed - i.e., the execution sequence 
is G400/F389/G400, performed twice (see Figure 25). The rotor hub 
load excitation vector, rotor impedance matrix, and HHC matrix and 
stress matrices are not formed. The check of the formation of 
these matrices with variable inflow obtained from the coupled 
G400/F389 program is obtained in the ARES model application 
described subsequently. 

Tai l Option Program (RIEVA) 

Case 1 verifies the calculation of harmonics of tail loads when 
rotor induced far wake data derive from a file to which inflow 
harmonics are transmitted by F389. Figure 32 shows the empennage 
airload variation with time. 

Case 2 is similar to Case 1 except that the effects of the far 
wake are not included. 

Rotor Aeroelastic Analysis E927 

Case 1 has the rotor attributes described under Base Program check 
Case 12. Rotor M, C, and K matrices were transmitted to the Base 
Program. 

Base Program Plot Package 

Case 1 contains file data transmitted to a plot file by the Base 
Program. The file data were calculated for the six degree-of- 
freedom isolator system described as Base Program check Case 11. 
Figures 15, 16, 19, and 20 illustrate plots that may be obtained 
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(see Reference 18 for plot variants). 


WABA T Geometry/Plot Program (APT/Plot) 

Case 1 checks the formation of panel geometry in the Automated 
Paneling Technique (APT/Plot) program from a hypothetical fuse- 
lage. The fuselage consists of three sections, each made up of 
several panels. The sections consist of a nose cone, center 
section cylinder, and reversed nose cone, which completes the 
fuselage. The APT/Plot program assembles the three sections into 
a body resembling a fuselage. The resulting geometry data are 
stored for input to the WABAT potential flow program or the plot 
module of APT/Plot. 

Case 2 employs stored geometry data to check the plot module of 
the APT/plot program. The plot file data contain panel geometry 
data for a typical fuselage. Plotted results are shown in Figures 
27 to 29. The results show the capability of the program to 
permit viewing of the geometry from any perspective, and the 
capability of the program to magnify regions of interest. 

Case 3 employs WABAT potential flow calculations to display flow 
directions (arrows in Figures 27 to 29) on the surfaces of fuse- 
lage panels. It was necessary to run the WABAT potential flow 
module prior to obtaining the fuselage flow plots. The case 
illustrates the surface flow plotting function of the APT/Plot 
package . 

WABAT Potential Flow Program 

Case 1 uses APT/Plot file data to obtain the potential flow 
solution for the fuselage shown in Figures 27 to 29. Harmonics of 
fuselage inflow at the rotor are calculated and can be transmitted 
to F389. 
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USER DOCUMENTATION 


The following section lists the docxments which must be used to 
apply the programs contained in the Coupled Rotor/Airframe Vibra- 
tion Analysis package. A brief description is given of the 
purpose and contents of the documents . 

1. Cassarino, S.and Sopher, R. : Coupled Rotor/Airframe Vi- 

bration Analysis Program Manual. Volume I. User’s and 
Programmer's Instructions NASA CR-165891, 1982. 

The manual contains user instructions for the Base Program des- 
cribing input preparation and output interpretation. Input/output 
files are described for transmitting data to External Programs. 
Program structure and subroutine functions are described. Appen- 
dixes describe the structure, file functions, and usage of the 
coupled G400/F389 program, and the E927 program. Descriptions are 
given of procedures for modifying the Base Program. Sample test 
cases are listed and described. 

2. Volume 2 of the above contains line printer input and output 
samples for the Base Program, and coupled G400/F389 and E927 
programs . 

3. Bielawa, Richard, L. : Aeroelastic Analysis for Helicopter 

Rotors with Blade Appended Pendulum Vibration Absorbers - 
Mathematical Derivations and Program User ' s Manual . NASA 
CR-165896,1982. 

The mathematical basis of modifications to the G400 analysis 
comprising teetering rotor and pendulum absorber modifications and 
the method used to calculate the rotor impedance are described. 
The complete mathematical basis of the analysis prior to modifica- 
tions is described in Reference 8 listed also as item 4 below. 
Program hierarchy and subroutine functions are described. Input 
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preparation is described. 


4. Bielawa, R. L. : Aeroelastic Analysis for Helicopter Rotor 

Blades with Time-Variable, Nonlinear Structural Twist and 
Multiple Structural Redundancy - Mathematical Derivation and 
Program User*s Manual. NASA CR-2638, 1976. 

5. Egolf, T. A., and Landgrebe, A. J. : A Prescribed Wake Rotor 

Inflow and Flow Field Prediction Analysis. NASA CR-165894, 
1982. 

This manual describes the mathematical basis of the rotor induced 
variable inflow analysis, F389, functions of program subroutines, 
and input preparation. The analysis is a single rotor version of 
the analysis listed below as item 6, which provides further 
details of the mathematical basis of the method. 

6. Landgrebe, A. J. , and Egolf, T. A.: Rotorcraft Wake Analy- 

sis for Prediction of Induced Velocities. USAAMRDL TR 
75-45, 1976. 

7. Cassarino, Sebastian, J. , and Mouzakis, Timoleon; Bifilar 

Analysis User's Manual - Volume II. NASA CR-159228, 1980. 

The manual describes input preparation, program structure, and 
functions of subroutines of the E927 program including a bifilar 
vibration absorber. Modifications were made to the bifilar 
version of E927 program for the present contract to eliminate the 
bifilar capability and to configure the E927 rotor as a sub- 
structure for use by the Base Program. Appendix B of item 1 above 
describes changes users must make to E927 input to employ the 
modified program and specifies the file for transmitting data to 
the Base Program. Sample inputs and outputs are provided in item 
1 above. The mathematical basis of rotor aeroelastic analysis 
E927 is described in Reference 6, listed as item 8 below. 
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Johnston, R. A., and Cassarino, S. J.; Helicopter Rotor 
Stability Analysis, USAAMRDL TR-75-40, U.S. Army, January 
1976. 


9. Gangwani, Santu T. : A Doublet Lattice Method for the 

Determination of Rotor Induced Empennage Vibration Airloads 

Analysis Description and Program Documentation. NASA 
CR-165893, 1982. 

The manual describes the basis, program structure, subroutine 
functions, and input preparation for use of the Rotor Induced 
Empennage Vibration Airloads (RIEVA) program. Applications are 
described and discussed. The file for accepting variable inflow 
data from the F389 program is specified in item 1 above. 

10. Studwell, R. E. : User’s Manual for the Coupled Rotor/Fuse- 

lage Vibration Analysis Graphics Package. NASA CR-165897, 
1982. 


The manual describes the user responses to plot package messages 
for plotting results transmitted by the Base Program. Plot option 
results are illustrated. Program function and subroutine struc- 
ture are described. The plot data transmission file format is 
specified in the text of item 1 above. 

11. Studwell, R. E.: User’s Manual for the Automated Paneling 

Technique (APT) and the Wing Body Aerodynamic Technique 
(WABAT) Programs. NASA CR-165895, 1982. 

This manual describes the preparation of input for the WABAT 
geometry package (APT/Plot) and the WABAT potential flow program. 
The manual describes the user responses to plot fuselage geometry 
and fuselage surface flows. Sample problems are illustrated. 
Input data requirements for the geometry package (APT/Plot) are 
described. The basis of the automated paneling technique is 
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described in Appendix A of Reference 13 (see item 12 below). The 
potential flow analysis basis is summarized in Appendix B of 
Reference 13 (see item 12 below). 

12. Sheehy, T. W. : A Method for Predicting Helicopter Hub Drag. 

(USAAMRDL Contract DAAJ02-74-C-0050 . ) USAAMRDL-TR-75-48 , 

1975. 


CORRELATION OF THEORY AND TEST 


In addition to development of the Coupled Rotor/Airframe Vibration 
Analysis, this contractual effort included a limited correlation 
study. The correlation study involved using wind tunnel data from 
a dynamically scaled UH-60A rotor system coupled with the NASA 
Langley aeroelastic rotor experimental system (ARES) powered test 
stand. Tests were conducted with and without higher harmonic 
control (HHC) employed. Model balance loads, rotor bending 
moments and vibrations in the ARES structure were measured. A 
shake test was performed to obtain the modal attributes of the 
system for input to the Coupled Rotor/Airframe Vibration Analysis. 

Shake Test of the ARES Model 

A shake test was conducted in the NASA Langley Transonic Dynamics 
Tunnel (TDT) to determine the support system modal properties 
required for vibration analysis predictions. The support system 
was provided by the Langley Research Center Aeroelastic Rotor 
Experimental System (ARES) (see Reference 19). The test procedure 
consisted of exciting the model by an electromagnetic shaker in 
several different directions. Random excitations were applied to 
the model. Data were recorded on 14 track FM magnetic tapes and 
consisted of outputs from eight hub accelerometers, three body 
accelerometers, six rotor balance channels, and the force load 
cell . 
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The shake test data were processed using a dual channel digital 
analyzer (HP 5420) to obtain transfer functions - e.g. accel- 
eration/applied force harmonics - of the measurements with respect 
to the applied forces . The ARES modal attributes presented in 
this section and those used in the correlation study are based on 
the random excitations. Table 8 lists the characteristics of nine 
ARES modes covering a frequency range from 21 to 73 Hz. The 
frequencies of interest are located near 45 Hz which is 4/rev. for 
the four bladed model rotor. The natural frequency, generalized 
mass, percent critical damping, six-component mode shape at the 
hub, and three translational components of modal response at the 
airframe accelerometer location for each mode are listed in Table 
8 . 

Modes 1, 2, 3, 7, 8 and 9 were obtained from a vertical off center 
excitation of the hub. A centered lateral excitation run defined 
the properties for modes 4 and 5, A centered longitudinal excita- 
tion run defined mode 6. The generalized masses for each of these 
nine modes were determined from single degree of freedom analyti- 
cal models of the system. The damping in each mode was obtained 
by peak fit using the HP 5420. Imaginary components of the 
transfer functions were used to define the mode shapes at the 
rotor hub and airframe . 

Wind Tunnel Test Program 

The test program was conducted in the NASA Langley Transonic 
Dynamic Tunnel (TDT). The model rotor was installed on the 
Langley Research Center Aeroelastic Rotor Experimental System 
(ARES), Reference 19. Figure 33 shows the ARES model with the 
test rotor. 

A four-bladed articulated rotor (1/6 scale) was installed on the 
ARES model. The dynamically scaled blades were based on UH-60A 
blades. The same blades were used for baseline studies in 
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Tunnel with Model Test Blades. 































Reference 20 which contains details of the blade properties. The 
blades were scaled for operation at full-scale Mach numbers in 
Freon. Model Reynolds numbers were 40 percent of full-scale 
values. The blades have a 56.22-inch radius, a 3.625 inch chord, 
and a -12 degree twist. The blades have an SC1095 airfoil section 
inboard of the 50 percent and outboard of the 85-percent radius 
stations. The SC1095R8 airfoil is used between the 50- and 
85-percent span stations. 

To achieve vibration reductions with HHC, various types of higher 
harmonic controllers were tested. The control algorithms are 
described in References 11 and 21. Comparisons of theoretical 
predictions obtained with the Coupled Rotor/Airframe Vibration 
Analysis were made with test data for the adaptive closed-loop 
controller listed as controller 2 in Table 1 of Reference 21. In 
this controller, the T matrix was identified and updated on-line 
and "cautious" terms were provided (see Reference 21). 

Application of the Analysis 

The Base Program was applied in its forced response solution mode, 
FRl . The forcing frequency used in this solution mode for the 
present application was 44.7 Hz which is the 4/rev. frequency for 
the four bladed rotor model at its operating speed of 670 RPM. 

Rotor Representation . The rotor effects were represented by 
substructure RES which models a multibladed aeroelastic rotor in 
forward flight or hover. This component reads the rotor hub 
excitation force and moment vector and the rotor impedance matrix 
and HHC and blade stress matrices computed by rotor aeroelastic 
analysis G400. In this application G400 was executed iteratively 
with External Program F389 to account for rotor induced variable 
inflow. The iterative sequence was G400-F389-G400 - i.e., one 
coupled cycle of these programs was run (see Figure 25). 


127 


The rotor was analytically trimmed to achieve the test values of 
lift and shaft angles for advance ratios from 0.20 to 0.40. Blade 
properties were taken from Reference 20. The resulting blade 
frequencies are listed in Table 9. Two flatwise, one edgewise, 
and one torsion mode(s) were employed in the computations. 
Full-scale Reynolds number data for the SC1095 and SC 1095R8 
airfoils were used in rotor aeroelastic analysis G400. The thrust 
coefficient C^/a was equal to 0.07 for all cases. 

Rotor Impedance Sensitivity to Perturbation Amplitudes . - Ex- 
ternal Program G400 obtains the rotor impedances with respect to 
hub motions and HHC angles. The perturbation amplitudes were 
selected to achieve accurate values of the numerically obtained 
impedance elements. Applications showed that there was a practi- 
cal limit on the size of the perturbation amplitudes for accurate 
results. Accurate results were verified by comparisons of G400 
numerical results with analytically derived results for rotor 
impedance and HHC control matrix elements . 

Figures 34 and 35 show sensitivity of the z^^(l,l) and h^^jj(3,3) 
elements of the rotor impedance and h^j matrices to perturbation 
amplitude. Based on these figures, the perturbation amplitudes 
selected were 152 m/sec^ (500 ft/sec^) for hub linear acceler- 
ation, 50 rad/sec2 for hub angular acceleration, and 0.10 degree 
for HHC angles . 

Support System Representation . - The ARES modal attributes listed 
in Table 8 were incorporated into the analysis by the use of modal 
substructure MSI. Each ARES mode was represented by one MSI 
substructure. The connection point between the rotor hub and 
support system was defined as the first node. The airframe 
location where vibrations were measured was taken as the second 
node. The airframe node required the use of a generalized force 
GFl type of substructure. 
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TABLE 9 


- CALCULATED BLADE NATURAL FREQUENCIES FOR 
THE MODEL ROTOR (PER REV) 


Mode No. 

1 

2 

3 

Flatwise 

2.7 

4.7 

7.7 

Edgewise 

5.0 

14.2 


Torsion 

6.0 

14.9 



TABLE 10 - COMPARISON OF ARES AIRFR?^ AND MODEL ROTOR 
BLADE NATURAL FREQUENCIES 


ARES SHAKE TEST DATA 
(REFERRED TO FIXED SYSTEM) 


CALCULATED 
BLADE MODES 

(REFERRED TO MOVING SYSTEM) 
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Specified Higher Harmonic Control » - Four values of pitch ampli- 

tude were specified to determine the effects of non-optimal 
control angle inputs on vibrations and loads. These four cases 
were obtained by specifying = + 1®, and -1°, and e^g = + I®, 

and -1®. These specified inputs are non-optimal because they are 
not determined from the requirement that they must minimize 
vibrations . 

Optimal Higher Harmonic Control . - The analytical optimal control 

state is derived from equation (93), which is essentially the same 

as the algorithm presented in reference 11. The control angle 

vector has six components listed in equations (78) and 79) - i.e. 

cosine and sine components of 3-, 4-, and 5-per-rev. inputs. 

Three airframe vibrations (longitudinal, lateral, and vertical) 

were selected for inclusion in the performance index, J, to be 

minimized (see equation (90)). Elements of the weighting matrix 

w were 20,000, 20,000, and 150,000 - i.e., vertical vibration 

z 

reduction v/as given more weight than rotor inplane vibration 
reductions. Diagonal elements of the control weight matrix used 
to limit the amplitudes of HHC control inputs were each equal to 
.01. Analytical control inputs were found to be less than 1.5® 
for all conditions except |j = 0.5 where 0^ was predicted equal to 
3®. The test HHC inputs were limited to a maximum of 1.5®. The 
analytical controller employed an analytical determination of the 
T matrix (see Higher Harmonic Control for Vibration Reduction). 
Noise was not represented. 

Correlations for Baseline Conditions 

Airframe Vibrations . - Comparison of predicted vibration levels for 
a three mode and nine mode solution versus ARES test data are 
shown in Figure 36. Both solutions adequately predict the trend 
of increasing vibration levels with increasing advance ratio. The 
modes used for the three mode solution are those modes with 
frequencies closest to the rotor 4/rev. which should be the primary 
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contributions to vibration levels . The nine mode solution was 
employed in an attempt to improve the correlation with test data 
and assess the sensitivity to modal inputs. As would be expected, 
the additional modes resulted in higher absolute values of vibra- 
tion levels, however, despite significant fuselage modal partici- 
pation of modes 2, 3, 8, and 9, the large contribution to vibra- 
tion levels was not expected. In fact, inclusion of the addi- 
tional six modes results in predicted vibration levels signifi- 
cantly higher than measured values . All modes used were obtained 
from random shake test excitations which historically have proved 
less accurate than swept sine results. Since the vibration 
predictions appear to be very sensitive to modal characteristics, 
a higher level of accuracy in representing the fixed system modes 
is desirable. 

Rotor Blade Moments . - Blade vibratory bending and torsion moment 
trends with airspeed predicted by the analysis (see Figure 37) 
were examined. Edgewise moments are overpredicted whereas the 
torsion moments are underpredicted by more than 50 percent. Blade 
moment variations with radius shown in figure 38 for 0.30 also 
have the same behavior. 

The torsion test data in Figure 38 appear to contain errors. Test 
data for C^/a = .075 and }j = 0.3 derived from Reference 20 are 
plotted in Figure 39. The test flatwise and edgewise moments show 
a slight increase in comparison with the curves in Figure 38. The 
torsion correlation is good, suggesting that the torsion test data 
in Figure 38 contain errors. 

Figures 40 and 41 show the harmonics of the blade moments at mid 
span and m= 0.30. Figure 41 contains test results from Refer- 
ence 20. The flatwise moment harmonics are predicted adequately 
by the analysis . 
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Figure 37. Variations of Blade 
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Figure 38. Radial Distributions of Blade 
Moments at 0.3 Advance Ratio. 
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The edgewise harmonics correlate satisfactorily except at 5- and 
6 -per- rev. The analysis predicts a large 5/rev, moment and a 
small 6/rev. moment. The large 5/rev. component accounts for the 
overprediction of the edgewise moments, (see Figures 37 and 39). 
This 5/rev. amplification is probably due to resonance of this 
mode with 4/rev. rotor hub excitation, as explained below. 

Pitch link load comparisons were similar to the torsional moment 
curves in Figure 37. 

The conclusion is that the blade moment trends are adequately 
predicted by the analysis- Absolute values appear to depend on 
how accurately the blade parameters are modelled, especially, 
blade frequencies - 

Hub/Balance Force Correlations . - Analytical 4-per-rev. hub force 
and balance test force trends with airspeed are shown in Figure 
42 . A direct comparison of hub loads is not made because the test 
data consist of balance loads while the analysis data consist of 
hxib loads. A dynamic calibration matrix is required to relate the 
4-per-rev. loads at the model hub to the 4-per-rev. loads sensed 
by the balance. The calibration matrix was not available in this 
correlation study. Figure 42 also shows that the longitudinal hub 
force from analysis has a greater sensitivity to airspeed than the 
longitudinal balance force. The lateral and vertical hub force 
trends with airspeed correlate fairly well with the corresponding 
balance loads . Since no direct comparisons of hub loads are 
available, no conclusions can be made on the accuracy of the 
analysis, except that the analysis appears to show the proper 
trend of hub loads with airspeed. 

Effect of Hub Displacements . - An area of interest is the sensi- 
tivity of blade moments and hub loads to support system flexi- 
bility. If the support system modal attributes are such that 
blade loads are substantially different from those with a fixed 
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Figure 42. Comparison of Variations with Advance Ratio of 
Predicted Hub Loads and Measured Balance Loads 
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hub, then the blade analysis should properly account for effects 
of hub displacements. 

Figure 43 shows that blade edgewise moments are affected by hub 
flexibility, and that hub motion increases the alternating edge- 
wise moment at mid span by 50 percent. Table 10 reveals the 
source of the increase in edgewise moments. The blade has an 
edgewise mode at 5/rev. which transforms into 4/rev. and 6/rev. in 
the fixed system. Only the 4/rev. component is transmitted to the 
fixed system for this 4-bladed rotor system. This transformed 
4/rev. edgewise mode is at resonance with the 4/rev. rotor excita- 
tion and may couple with the fixed system lateral roll mode at 
4.4/rev. This speculation is reinforced by Figure 44 which shows 
that the in-plane hub forces are more sensitive to hub displace- 
ment than the vertical hub force. 

The deduction is that hub participation can have a substantial 
effect on blade moments and hub loads. The accurate prediction of 
hub loads is sensitive to accurate characterization of support 
system modal attributes and rotor blade frequencies. The lack of 
agreement of edgewise moments with test data is likely to be due 
to inaccurate characterization or rotor blade natural frequencies 
in the analysis. Because rotor loads are the primary excitation 
sources of coupled rotor/airframe systems, these sensitivities 
show the importance of accurate modal characteristization of both 
the rotor and airframe for vibration prediction. 

Correlation with HHC Non-optimal Inputs 

The Coupled Rotor/Airframe Vibration analysis was executed to 
predict airframe vibrations, rotor bending moments and hub loads 
utilizing arbitrary higher harmonic control inputs. Comparison 
with wind tunnel measurements for the same inputs is discussed 
below. 
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Airframe Vibrations . - Figure 45 shows the effect of 3/rev. HHC 
pitch inputs on airframe 4/rev. vibration. Also shown is the 
analytical estimate of vibration obtained from the Coupled Rotor/ 
Airframe Analysis with HHC and a three mode representation. 
Absolute values are underpredicted and the trend in vibration 
reduction with HHC application is not predicted uniformly. The 
nine mode solution results shown in Figure 46 in general agree 
better with the test data and shows uniformly higher vibration 
levels than the three mode solution. 

Figure 47 shows vibration vector corrrelations for the airframe 
4/rev. longitudinal vibrations. For all four HHC input values, 
the analysis is able to predict the phase changes due to HHC 
inputs . Amplitude comparisons are the same as in Figure 45 . 

Rotor Blade Moments . - Figure 48 shows the effect of HHC inputs on 
blade bending moments. This figure indicates that bending 
moments, especially edgewise moments, are sensitive to the HHC 
inputs. Harmonic analysis of the moments (see Figure 49) indi- 
cates the differences are primarily due to large 3/rev. flatwise 
and 6/rev. edgewise contributions for the HHC applications. The 
predicted blade moment distribution for baseline and HHC input 
case are provided in Figure 50. The sensitivity to HHC inputs 
observed in the measured data is not predicted. 

The sensitivity of the measured blade loads to a 3/rev. input may 
be due to impurities in the 3/rev. pitch input and the proximity 
of the test blade edgewise mode to 6/rev. The impurities are 
assumed to give rise to a 6/rev. excitation which excites the 
6/rev. edgewise mode. This observation was made in Reference 11. 
The analysis which models a 3/rev. HHC input as a 3/rev. pitch 
change is not able to predict this effect. 
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Figure 45. Effects of Specified HHC on Amplitudes of 
Airframe Vibrations - Three Mode Solution. 
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Figure 46. Effects of Specified HHC on Amplitudes of 
Airframe Vibrations - Nine ,Mode Solution. 
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Figure 48. Effects of Specified HHC on Measured 
Blade Moment Distributions. 
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Figure 49. Effects of Specified HHC on Harmonics 
of Measured Blade Moments. 
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TORSION MOMENT 




Hixb/Balance Force Correlations ♦ - Figure 51 shows side-by-side 
comparisons of predicted hub forces with test balance forces. 
Variations with HHC inputs are not well correlated. Because no 
direct hub load comparisons are made, no conclusions are drawn. 

Optimal HHC Solutions 

Analytical and test results are compared for the effects of 
optimal values of HHC inputs on airframe vibrations, blade 
moments , and hub loads . 

Airframe Vibration . - Predicted and measured airframe vibrations 
for the baseline model and optimal HHC inputs as functions of 
advance ratio are shown in Figure 52. The analysis predicts the 
proper trends with HHC application and advance ratio increases for 
both the vertical and longitudinal accelerations. The test data 
indicate that the HHC inputs are not effective in reducing lateral 
vibrations in contrast to the analytical predictions. 

Blade Bending Moments . - Both analysis and test show that for the 
conditions selected, the spanwise blade bending moment distribu- 
tions are not sensitive to HHC inputs. Figure 53. The disagree- 
ment in the torsion moments is probably due to errors in the test 
data as noted in the correlation of baseline conditions. 

Figure 54 shows that predicted flatwise and edgewise moments are 
consistent with measured values and that the trends with advance 
ratio are predicted correctly. The trend for torsional moments is 
consistent but the predicted magnitudes are consistently lower 
than measured. 

Figure 55 shows comparisons between test results and predicted 
blade moment harmonics at midspan and 0.30. The torsion 
moments will not be discussed because the test data probably 
contain errors. The flatwise and edgewise harmonics of moments 
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Figure 56. Effects of Optimal HHC on Comparisons of 

Predicted Hub Loads and Measured Balance Loads 




correlate fairly well. The overprediction of the 5/rev. and 
underprediction of the 6/rev. edgewise components may be due to 
the input blade edgewise mode frequency being 5/rev. while the 
actual frequency is closer to 6/rev. , as discussed in the pre- 
ceding text. 

Hub/Balance Force Comparisons . - Figure 56 shows that the analysis 
predicts a large effect of the HHC inputs on the longitudinal hub 
force which is not evidenced in the measured balance loads. The 
effects on lateral and vertical forces are more consistent but the 
measured balance vertical forces are more sensitive to HHC inputs 
than the predicted effect on vertical hub forces. 

CONCLUDING REMARKS 


1. A Coupled Rotor/Airframe Vibration Analysis has been 
developed to serve as a preliminary design tool for pre- 
dicting helicopter vibrations and a research tool to 
quantify the effects of structural properties, aerodynamic 
interactions and vibration reduction devices on vehicle 
vibration levels. 

2 . The analysis developed employs a substructure method which 
provides an efficient stand-alone mode operation of a Base 
Program while retaining a high level of aerodynamic and 
dynamic sophistication by linking output of External 
Programs . 

3 . The Base Program operating in a stand-alone mode can 
efficiently be used to study design changes in fuselage 
properties and vibration reduction devices . 
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4. 


The substructure basis of the Base Program is well suited to 
allowing the Base Program to introduce data from user 
provided External Programs of the rotor system. The sub- 
structure basis also will facilitate growth of the Base 
Program to include additional internal components. 

5. The Rotor/Airframe coupling in the analysis is accomplished 
using an impedance matching technique. 

6 . The current External Programs supporting the Base Program 
provide the capability to model articulated, hingeless and 
teetering rotor systems, higher haimionic control inputs, 
main rotor wake/empennage excitation sources and hover, 
transition and forward flight regimes. 

7. A limited correlation study with a dynamically scaled wind 
tunnel model indicates: 

a. The analysis predicts the proper trends in vibration 
levels as a function of advance ratio. 

b. The analysis predicts the proper trends in vibration 
levels and blade bending moments with application of 
higher harmonic control . 

c. Absolute values of predicted vibration levels and blade 
bending moments are sensitive to airframe modal char- 
acteristics and were not uniformly consistent with 
measured values . 
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APPENDIXES 


Appendix A - Illustration of System Synthesis 


Substructure Terminology 


Connection 
node coordinate 


Connection 
node coordinate 



Unconnected Coordinate 


Element 1 


Connection Node 1 


Connection Node 1 


Element 2 


Sketch (Al) 


Substructure Equations of Motion 


Element 1 
Coordinates 

Properties 



(Al) 


(A2) 


(A3) 


160 




Ex1:ernal force 
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[O] 
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force 




(A14) 


Reaction 
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Assembled equation employing dependent coordinates 


(A16) 


Mapping From Dependent Coordinates To Independent Coordinates 


Component displacements are equated at connection node 1 
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System Synthesis 


Properties of assembled system referred to independent coordinates 


M = 

K = 

F = 

■ 1 -. 

because 


m^ O 
O O 



(A22) 


(A23) 


(A25) 


Assembled system equation of motion 



(A27) 


MXj + K Xj = O (A28) 

Recoverx of Component Responses and Reactions 
Components responses 

Xjj = pXj (A29) 

Xp = pXj (A30) 

Xp = pXj. (A31) 
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(A33) 


Component reactions 

Rd = - ‘"d'Sd ^ 
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Rotation Transformation 
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For no rotation B = I 
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Organization of Coordinates and X^^ 


Elem 4 
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Composite Transformation 


For the case where substructure and global axis orientations are 
the s ame 
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Appendix B - Substructure Equations 
Fixed Absorber FAl 

See Figure 5 for schematic of absorber. 
Coordinates 


No. 

X 

1 

6 

2 

u 

3 

V 

4 

w 

5 

8l 

6 


7 

«3 


Matrices for the fixed absorber are shown in Figure B1 . 

The force vector is (see Figures 2 and 5) 

{f} = {O} (Bl) 


The reaction load vector (see Figure 2 for positive load direc- 
tions) is 


{r} ^ = {O, r^^. 

Properties are 


^yl' ^zl' ^mxl' ^myl' ^mzl^ 


(B2) 


m - mass 


I , I , I - principal mass moments of inertia of absorber 
4 - damper constant, percent critical 

c = 2tuj^m, damping coefficient (B3) 

k - stiffness 

u>^ = (k/m)^, natural frequency of absorber (B4) 
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Generalized Force GFl 


See Figure 6 for schematic of svihstructure . 
Coordinates 


No. 

X 

1 

u 

2 

V 

3 

w 

4 

®1 

5 

CD 

6 

83 


The matrices for the GFl substructure are 


[m] 

= 

[O] 

(B5) 

[c] 

= 

[ 0 ] 

(B 6 ) 

[k] 

= 

[O] 

(B7) 


The specified load vector at the connection node is 


~ ^^xl' ^yl' ^zl' ”^xl' ”^yl' "^zl^ 

The reaction load at the connection node is 


(B 8 ) 


{r} - r^^, ^myl' ^mzl^ 


(B9) 
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Beam Model BMl 


See Figure 7 for schematic of beam. 
Coordinates 


' '1 

No. 

X 

No. 

X 

1 


7 



1 


2 

2 

V- 

8 



1 


2 

3 

w-1 

9 



1 


2 

4 

<®1>1 

10 

(®1>2 

5 

<®2>1 

11 

«®2>2 

6 

<®3>1 

12 

<«3>2 


Figure B2 illustrates the mass matrix of the system. Figure B3 
illustrates the stiffness matrix. The applied force vector is 


{f} = {0} 


(BIO) 


The reaction load vector 


^^^xl^l' ^^yl^l' ^^zl^l' ^^mzl^l' ^^myl^l' ^^mzl^l. 


^^zl^2' ^^mxl^2' ^^myl^2' ^^mzl^2 ^ 


(Bll) 


Directions of force components at a connection node are shown in 
Figure 2 . 
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Figure B3. Stiffness Matrix of Uniform Beam Segment^ 


Properties of the beam are 


m 

A 

E 


L 


YA 


ZA 


GJ 


El, 


El, 


mass 

cross-sectional area 
Elastic modulus 
Length 

principal area moment of inertia about y axis (see 
Figure 7 ) 

principal area moment of inertia about z axis (see 
Figure 7 ) 

torsional rigidity 

ElyA (B12) 

(B13) 


Constraint Substructure CNl 


Coordinates 

{q} ^ q2/- • • • q^jj} (B14) 

1 < M < 6 

Matrices 


[m] = 

[O] 

M X M 

(B15) 

[c] = 

[O] 

M X M 

(B16) 

[k] = 

[O] 

M X M 

(B17) 
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{f} 


{ 0 } 


M X 1 


(B18) 


Constraint substructure CNl is a special type of modal substruc- 
ture. Equation (20) is the basis of the constraint procedure. As 
an example, if it is desired to constrain - i.e. eliminate as 
coordinates - all physical displacements at a connection node 
except two - say, and 0^ - equation (20) takes the special form 


{X} 


Cy] . {q} 


1 

0 

0 



1 

0 


rA 

o 

0 


1 

H 

CD 

0 

0 


hA 

0 

1 






\ 0 ^ / 

o 

o 


\ 3 1 



— 


GLOBAL 

ORIENTATION 

AXIS 



(B19) 


This is obtained by assuming [E] = I and replacing elements of the 
modal matrix by zeroes where it is desired to constrain coord- 
inates, and ones where it is desired to retain coordinates. The 
modal coordinates q^ and are identical to v^, Q^. Analogous 
transformations apply to other selection combinations. 

Aeroelastic Rotor Type 2 (RE2) 

Rotor coordinates are listed in Table 4. Equations of motion are 
defined in reference 6. 
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Modal Substructure MSI 


Coordinates 


No. 

X 

1 

<^1 


Matrices 

[m] 

= [m^i] 

(B20) 


[c] 

= [c^i] 

(B21) 


[k] 

= 

(B22) 


{f} 

II 

0 

(B23) 


At a connection node, the generalized reaction load vector is 

T 

{Q}j^ = iy] {r} (B24) 


The modal vector at a connection node is 

{y} ~ ^^^xl^l^ ^^yl^l^ ^^yl^l^ ^^01^1' ^^02^1^^ ^^03^1^ (®25) 

The elements of the modal vector are displacements which have the 
positive directions shown in Figure 2. 

The reaction load vector at a connection node is 

{r} = ^yl' ^zl' ^mxl' ^myl' ^mzl^ 

The displacement vector at a connection node is 

{X} = (y) {q} 

T 

{x} = { u^, v^, w^, 0 ^, 0 ^/ 63 } local orientation 

axes 

If a physical substructure is not connected to a node of the modal 
substructure, obtain the physical displacements of the node by 
connecting a GFl substructure to the node with zero applied loads 
specified in GFl. 


(B26) 

(B27) 

(B28) 
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Properties 


m 

c 


11 

11 


k 


11 


4 


U) 


n 


generalized masss of tree normal mode 
generalized damping of mode 

generalized stiffness of mode 
damping of normal mode, percent critical 
natural frequency of the mode 
H 


(k^l/n.11 


) 


Horizontal Linear Bifilar BFl 


(B29) 


(B30) 


Coordinates 


No. 

X 

1 


2 

'‘b 

3 

u 

4 

V 

5 

w 

6 


7 

®2 

8 

83 


The multi -blade coordinates y and y are related to the 

c s 

bifilar arm lead physical displacements by 


y = (y^ cos»|j + y^ siniji)(2/N) (B31) 

See Figures 8 and 2 for positive directions of y and connection 
node displacements. Figure B4 shows the non-zero submatrices of 
the mass, damping and stiffness matrices. 
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The applied load vector is 


if] = {O} 


The reaction load vector is 


(B32) 


Figure 
loads , 


= { 
2 


^xl' ^yl' 
illustrates the 


positive 


^myl' ^mzl^ 
directions of 


the 


(B33) 

reaction 


Properties 


R 

r 

n 

Q 

N 


bifilar arm mass 

equivalent hinge offset 

equivalent pendulum length 

tuning ratio 

(R/r)h 

damping, percent critical 
rotation speed 

number of bifilar arms (N > 2) 


(B34) 
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Vertical Linear Bi filar BF2 


Coordinates 


No. 

X 

1 

0 

CQ. 

2 

Pc 

3 

Ps 

4 

U 

5 

V 

6 

w 
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8 

®2 

9 

®3 


The multi-blade coordinate bifilar arm flap displacements 
are related to the bifilar arm physical flap displacements by 

p = p /N + O cosij) + p sin4>)(2/N) (B35) 

o c s 

See Figure 9 for the positive directions of p and connection node 
displacements. Figures B5 to B7 show the mass, damping, and 
stiffness matrices. 


The applied load vector is 


{f} = {O} 


(B36) 


The reaction load vector is 


{r}*^ - {O, O, O, r^j^, r^^, r^^, ^myl' ^mzl^ 

Figure 2 illustrates the positive directions of the reaction 
loads . 


(B37) 
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Properties 


™b 

R 

r 

U) 


bi filar arm mass 
equivalent hinge offset 
equivalent pendulum length 
collective mode flap frequency 

ft (r+R)Vr 


4 - damping, percent critical 

ft - rotation speed 

N - number of bifilar arms (N > 2 ) 

Nodal Isolator ISl 


Coordinates 


No . 

X 

No- 

X 

1 


13 

^BF 

2 

Yt 

14 

^BF 

3 


15 

^BF 
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16 

CD 
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17 

<»2>2 

6 

"‘b 

18 

<®3>2 

7 

^BT 



8 

^BT 



9 

^BT 



10 




11 




12 

<«3)l 




(B38) 
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Figures B8 and B9 illustrate the non-zero elements of the 
mass, damping, and stiffness matrices for the isolator. The 
applied force vector is 


{f} = {O} 

The reaction load vector is 


{r} 


(B39) 


{ 0, 0, 

o 

o 

o 

o 







<W>1' 

^^xl^2' ^^yl^2' 

<^zl>2' 



^^mzl^2^ 




myl ^1 
mxl ^2 


(B40) 


Directions of force components at a connection node are shown in 
Figure 11 . 


Figure 11 illustrates the isolator with a single antiresonant bar 
weight, and Figure 12 is a schematic showing the isolator with the 
second antiresonant bar weight. 


Properties of the isolator are A, B, C, D - bar dimensions denoted 
by , B_ , C- , and , respectively, in Figure 12. 

Li Li Li Li 


m 


lx 


m 


ly 


mass of first isolator bar for displacement 

w^/g (B41) 

effective mass of first isolator bar for y^ dis- 
placement 

h 

(Wj^/g)/[(l- U)|y)2 + <2£iy«*3^y)2] (B42) 
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(D/C) 
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Figure B8. Non-zero Mass Submatrices for Isolator ISI. 







m 


Iz 


m 


2x 


effective mass of first isolator bar for z^ dis- 
placement 

(w^/g)[(l- m2^)2+(2£j^^U) (B43) 

effective mass of second isolator bar for x„ 

a 

di spl acement 

H 

(W2/g)/[(l- 


w^, “ antiresonant bar weights 


U) 


ly' “iz' ”2x 


to. 


ratios of natural frequencies of antiresonant 
bars to forcing frequencies for z^, and x^ 

displacements, respectively. 


^ly' ^Iz' 


2x 


damping ratios of antiresonant bars for y^, z^, 
and x^ displacements, respectively. 

l5 


lY 


effective mass moment of inertia of first anti- 
resonant bar about y axis. 


IZ 


1 H 

effective mass moment of inertia of first anti- 
resonant bar about z axis 




(B46) 
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I 


zY 


effective mass moment of inertia of second anti- 
resonant bar about y axis. 






(B47) 


1 


lY' 



mass moments of inertia of isolator bars. 


Stiffnesses listed below are illustrated in Figure 11 where they 
are prefixed with T. 


*^EY' ^EZ 


^X' ^TY' ^T2 


equivalent beam stif fnessesses for anti- 
resonant bar displacements x_- x„, y,.- y„, 

i5 /J. 13 J; 

z„-z„, respectively. 
a j. 

stiffnesses of springs between bar pivot 1 
and connection node 1 . 


^AX' ^AY" ^AZ 


stiffnesses of springs between bar pivot 2 
and connection node 2 . 


^0B' ^0T' %T 


equivalent beam stiffnesses for antiresonant 
bar rotations 0^, 0^, and respectively. 




<*T 

- X3>/C 

®T 



- 

®T 

= 

<^B 

- Yy)/A 


(B48) 

(B49) 

(B50) 
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Damper elements are in parallel with stiffness elements shown in 
Figures 11 and 12. Damper elements corresponding to stiffness 
elements listed above are obtained by replacing K by C. Damper 
values are 



^TX 

^^TX ^TX^*"tx 


(B51) 


Cty 

2^Ty ^tY^^^TY 


(B52) 


Sz 

^^TZ ^TZ^^^TZ 


(B53) 


^AX 

^^AX ^AX"^*^AX 


(B54) 


Sy 

^^AY^AY'^^AY 


(B55) 


Sz 

^Sz Sz"^^AZ 


(B56) 


Sb 

^^0B Sb^‘^0B 


(B57) 


St 

2^qT St^‘^0T 


(B58) 


^0T 

^Sb Sb^*"0B 


(B59) 

where 

etc . 

are damping ratios , 

and ^TX' 

etc . are 

frequency constants 

employed to define damper values , 


The isolator 

with 

two antiresonant bar 

weights (w^ and w ) is 

referred to as the 

3-D isolator, and the 

isolator with 

one anti- 

resonant bar 

weight (w_) is referred to as the 2-D 

isolator . 

Properties of 

the 2 

-D isolator are obtained from the 3- 

D isolator 


properties by equating to zero the following parameters : 

^2' ^2Y’ ^0B' Sb' St* 

The following parameters are equated to one: 

C(= Cj^), D(= Dj^). 
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Aeroelastic Rotor Type RE3 

See Assembly of Rotor Impedance Matrix for coordinates and 
equations of motion. 
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